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Abstract 

In recent years, a rapidly growing literature has focussed on the construction of wavelet systems to 
analyze functions defined on the sphere. Our purpose in this paper is to generalize these constructions 
to situations where sections of line bundles, rather than ordinary scalar-valued functions, are con- 
sidered. In particular, we propose needlet-type spm wavelets as an extension of the needlet approach 
recently introduced by 37. 38 and then considered for more general manifolds by [181 1191 [20] . We 
discuss localization properties in the real and harmonic domains, and investigate stochastic properties 
for the analysis of spin random fields. Our results are strongly motivated by cosmological applica- 
tions, in particular in connection to the analysis of Cosmic Microwave Background polarization data. 
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1 Introduction 

In recent years a rapidly growing literature has focussed on the construction of wavelet systems on the 
sphere, see for instance [3] , [51] [H [2] and the references therein. These attempts have been motivated by 
strong interest from the applied sciences, for instance in the areas of Geophysics, Medical Imaging and 
especially Cosmology/Astrophysics. 

As far as the latter are concerned, special emphasis has been devoted to wavelet techniques for the 
statistical study of the Cosmic Microwave Background (CMB) radiation data. These data can be viewed 
as providing observations on the Universe in the immediate sequel to the Big Bang, and as such they 
have been the object of immense theoretical and applied interest over the last decade [H]. In particular, 
experiments such as the WMAP satellite from NASA have provided high resolution observations on 
the "temperature" (i.e. intensity) of CMB radiation from all directions of the full-sky ([57]). These 
observations have allowed precise estimates of several parameters of the greatest interest for Cosmology 
and Theoretical Physics. Spherical wavelets have found very extensive applications here, see for instance 
([23l[9l[46l[8l|35l[36l|52]) and many others. The rationale for such a widespread interest can be explained 
as follows: CMB models are best analyzed in the frequency domain, where the behavior at different 
multipoles can be investigated separately; on the other hand, partial sky coverage and other missing 
observations make the evaluation of exact spherical harmonic transforms troublesome. The combination 
of these two features makes the time-frequency localization properties of wavelets most valuable. 
Besides providing measurements on the radiation intensity, experiments such as WMAP have also pro- 
vided some preliminary observations of a much more elusive physical entity, the so-called polarization 
of the background radiation. Polarization is a property of electromagnetic radiation whose physical sig- 
nificance is described for instance in ( [28[ [53|). see below for more mathematical discussion. So far, 
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empirical analysis of polarization has been somewhat limited, because this signal is currently measured 
with great difficulty. The situation with respect to CMB polarization data, however, will significantly 
improve over the next years, for instance by means of the ESA satellite mission Planck (expected to be 
launched in Spring 2009) which will take full-sky measurements of the polarized CMB sky (fST) with 
much greater precision. Moreover, both ESA and NASA are planning high sensitivity full-sky satellite 
borne experiments within the next f 0-20 years. 

Polarization measurements are of extreme interest to physicists for several reasons. Indeed, not only 
do they allow improved precision for estimates of physical parameters which are already the focus of 
CMB temperature data, but they also open entirely new areas of research. Just to mention a striking 
example, at large scales the polarization signal is expected to be dominated by a component related 
to a background of gravitational waves which originated in the Big Bang dynamics (in the so-called 
inflationary scenario, see ([12]). Detection of this signal would be an outstanding empirical validation 
of many Theoretical Cosmology claims, directly related to Big Bang dynamic models but deeply rooted 
in General Relativity. Taking into account the huge amount of polarization data which will be available 
in the next 1-2 decades, as well as the important cosmological information contained in these data, it is 
clear that suitable mathematical tools for data analysis are in high demand. While a large amount of 
mathematical statistics techniques have been developed for analyzing CMB temperature data, far fewer 
mathematical tools are available to analyze polarization data. 

From the mathematical point of view, as we shall detail below, polarization can be viewed as a (random) 
section of a line bundle on the sphere. Our purpose here is then to extend spherical wavelet constructions 
to the case where sections of line bundles, rather than ordinary functions, are considered. To the best of 
our knowledge, this is the first attempt to introduce wavelet techniques for the case where one deals with 
more general mathematical structures than ordinary (scalar- valued) functions on manifolds. We believe 
this area can be expanded into disciplines other than cosmology - for instance, tensor-valued random 
fields emerge naturally in brain imaging data ( [43j ) . 

In particular, in this paper we shall extend to the line bundle case the so-called needlet approach to 
spherical wavelets. Spherical needlets were recently introduced by [37l[33, and further developed, and 
extended to general smooth compact Riemannian manifolds in [181 119( I20j . In a random fields environ- 
ment, needlets were investigated by ([HH]), with a view to applications to the statistical analysis of CMB 
data; applications in the physical literature include [41], [33], [25], see also ([29], [TT] [M] [15] ) and [30l[34] . 
More precisely, we will focus on needlet-type spin wavelets, and we will argue below that they enjoy 
properties analogous to those of the usual needlets in the standard scalar case. In particular, we shall 
show below that needlet-type spin wavelets enjoy both the localization and the uncorrelation properties 
that make scalar needlets a powerful tool for the analysis of scalar-valued spherical random fields. More 
details on the plan and significance of our paper are given in the next section. 

The companion article [17], intended primarily for physicists, goes into further detail about the statistical 
uses of spin wavelets. One should be aware, however, that there are some notational differences between 
that article and this one. 

2 Plan of the Paper and Significance of the Results 

The spin s concept that we will use in this paper originates in the fundamental work of Newman and 
Penrose [39j . Writing for physicists, they say that a quantity rj defined on the sphere has spin weight 
s S Z, provided that, whenever a tangent vector m at a point transforms under coordinate change by 
m' = e*'^m, the quantity i], at that point, transforms by rj' = e'^^^rj. They then develop a theory of spin 
s spherical harmonics, upon which we shall build to produce spin wavelets. 

In section 3, we put the notion of spin s into acceptable mathematical language. We do this in an 
elementary manner which remains very close to the spirit of the work of Newman and Penrose. For 
another approach, using much more machinery, see [13j . 

Let N be the north pole (0,0,1), let S be the south pole, (0,0,-1), and let Uj be 5^ \ {N, S}. If 
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R e S'0(3), we define Uj^ = RUi. On [// we use standard spherical coordinates {9,4>) < 9 < it, 
— TT < </> < tt), and analogously, on any Ur we use coordinates {6B.,(t'R) obtained by rotation of the 
coordinate system on Ui. At each point p of Ur wc let pji{p) be the unit tangent vector at p which is 
tangent to the circle 9r = constant, pointing in the direction of increasing (^j^. (Thus, for any p G [//, 
Pnk^v) = ^*v\Pikp)\) If P e f/i?i n Ur^, we let V'pfla.Ri be the angle from Pr^(j)) to Pr^{v). 
Now say C 5^ is open. We say that / ~ {!r)r.^so{z) G C'^(^)i or that / is a smooth spin s function 
on rj, provided aU ]r e C°°{Ur n fl), and that for aU i?2 G 50(3) and all p G L/fli n [/^^ n fl, 

fRM^e^^^-fRAp), (1) 

where V' — i'pR2Ri- If; ^^Yi -^i = ^2 = then heuristically //j is // "looked at after the coordinates 
have been rotated by i?" ; at p, it has been multiplied by e^^^ . The angle -0 would be clearly be the same 
if we had instead chosen Pr(jp) to point in the direction of increasing 9r, say, so this is an appropriate 
way to make Newman and Penrose's concept precise. Note that if s = 0, we can clearly identify C^(r2) 
with C°°(17). 

For any s, we can identify C^(ri) with the sections over Vl of the complex line bundle obtained by using 
the e*''^r'^2«i as transition functions from the chart Ur^ to the chart Ur^; then Jr is the trivialization of 
the section over Ur fl VL. This is the point of view that we take in section 3, since, as is often the case 
in mathematics, certain properties are clearer, and more easily verified, if one uses the coordinate-free 
line bundle point of view. However, right now, for the reader's convenience, we state our results without 
reference to line bundles. Line bundles are rarely mentioned explicitly in this article after section 3. 
We may define Ll{n) (resp. Cs{^)), by requiring that the Jr be in L'^{Ur{^^1) (resp. C{Ur{^^1)) instead 
of C°°{Ur n fi). There is a well-defined inner product on L^(il), given by (/, = {fR,gR); clearly this 
definition is independent of choice of R. A key observation of section 3 is that there is a unitary action of 
50(3) on L1{S'^), given by / ^ which is determined by the equation {f^)i{p) ~ fR{Rp). We think 
of as a "rotate" of /. 

In section 4, we explain the spin s theory of Newman and Penrose, within the rigorous framework we 
have just outlined. Most of our arguments are very close to theirs. 
For / smooth as above, following Newman and Penrose we define 

dsnfR = -(sin..)^ + -^^) (sin..)-V., 

and we show that the "spin-raising" operator d : Cf'{Vt) C^^IQ) given by {df)R = Osr/r is well- 
defined. We also define Osr/r by OsrJr = QsrJr, which leads to the spin-lowering operator d : (7^(17) — > 
0^]^(il) given by {df jR — Osr/r- We show that d and d commute with the actions of 50(3) on smooth 
spin functions. 

Now, for I > |s|, let hs = [{1 + s)\/{l - s)\]^^^ for s > 0, h, = [{I - 3)1/(1 + s)\]^/^ for s < 0. Let 
{Yim : I > 0, —I < m < 1} he the standard basis of spherical harmonics on S^. Following Newman 
and Penrose, for I > \s\, we define the spin s spherical harmonics by sYi,„ = d'^Yim/bis for s > 0, 
sYim — {—dyYim/bis for s < 0, and we show that {si^m : I > Nl,— ^ < to < /} forms an orthonormal 
basis for L^(5^). There is a (relatively) simple explicit expression for the sYimi', it has the form 

sYi,ni{9,cf>)^ syim{9)e'"'^ (2) 

for a suitable function sUim- (See (|60|l for the explicit formula.) li R I, gyimR does not have a simple 
expression; but sYi^j^{9r, equals syim[9 R)e'^'^'^ in the coordinates on Ur. (Here again sl^^ denotes 
the "rotate" oi gYim.-) 

There is also an analogue As of the spherical Laplacian for spin s functions; specifically, if s > 0, we let 
As — —(9(9, while if s < 0, we let As = —dd. Then Aq is the usual spherical Laplacian. For I > |s|, 
let His =< sYira : -I < m < I >. Also, for I > \s\, let Xu = {I - s){l + s + 1) if s > 0, and let 
Ajs = {I + s)(' — s + 1) if s < 0. Then His is the subspace of 0^(5^) which consists of eigenfunctions of 
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for the eigenvalue Xis- This of course generalizes the situation in which s = 0. All of the aforemen- 
tioned facts follow by making the arguments of Newman and Penrose rigorous. 

In sections 5, 6, 7 and 8, we present new results. In section 5, we define operators with smooth kernels, 
from C^{S^) to itself. We say K — (-ftrfl',_R)fl',fleso(3) is such a smooth kernel if each 
Kr',r eC°^{S^xS^), and if X^^^^, (x, = e''<^^'-^^KR,^nix,y) for all x € Ur^'HUr' and y G Ur.DUr, 
where = ^P^r/ r, and ip = ipyRR-^ ■ It is then evident that we may consistently define an operator 
/C : Cr(52) ^ C^iS^) by (/C/)i^'(x) = J^, KR,,R{x,y)fR{y)dSiy), for aU x € Ur,; we cah /C the 
operator with kernel K. (These operators may be identified with operators with smooth kernels acting 
on sections of the aforementioned line bundle.) 

Next we show that there is a good notion of spin s zonal harmonic. To see how this arises, we first show 
(in section 3) that, if / S (7^(5'^), then the function e^^'^fi{9, </>) extends continuously from Uj to Uj U N, 
and the function e~^'''^ fi{9, (p) extends continuously from Uj to Ui U S. (This is shown by considering 
Jr for those R with N, S S Ur.) Thus there is a well defined hnear functional L : Cs{S^) C, given by 
Lf — lime^Q+ e^^'^ fi{9^4>). It is evident from ^ that L{sYim) must be zero if m ^ ~s. This indicates 
that sYi,-s could have a special status; in fact, it is a multiple of the spin s zonal harmonic, in the sense 
of the following results, which we prove in section 5. 



Liemmas on Zonal Harmonics De/ine .Z, = (-1)^ ^/{2l + I)/ An sYi-s- Then: 

(a) L{sZi) = {21 + l)/47r. 

(b) Foranyfenis,Lf={f, ,Zi). 

Next, let Pis be the projection in L?{S'^) onto Tiis, so that Pis has kernel K^^ , where for any rotations 
R,R', K'jl, R{x,y) = Y.''^^_i sYimR'ix)sYimR{y) for x € Ur' , y £ Ur. Then: 

(c) IfxeUiD Ur, and N e Ur, then sZii{x) = e*^('^i-'''2)i^]|,_^(x, N), where V-i = ^xRfi, and ip^ «s 
the angle from Pr^N) to d/dy at N. 

(d) Ifxe Ur, then K^^ j^{x,x) = {21 + l)/47r. 

In section 6, we begin our study of spin wavelets. Before explaining our construction, we need to explain 
the ideas which are used in the spin case on the sphere, and on more general manifolds. The references 
here are [37], [38], [18], [H], [20]. 

A word about notation: in sections 3-5 we often use the variable "/" to denote a spin function on the 
sphere. In sections 6-8 the variable "/" will be reserved for another purpose. 

Specifically, say / S iS(R"'"), / ^ 0, and /(O) = 0. Let (M,!?) be a smooth compact oriented Riemannian 
manifold, and let A be the Laplace-Beltrami operator on M (for instance, the spherical Laplacian if M 
is 5^). Let Kt be the kernel of f{t^A). Then, as we shall explain, the functions 

wt.x{y) ^ Kt{x,y), (3) 

if multiplied by appropriate weights, can be used as wavelets on M. In case M = 5" and / has compact 
support away from the origin, we shall say these wavelets are "needlet-type" ; the theory of needlets was 
developed in [37], [38]. (Actually, the definition of needlet in [37], [38] is slightly different from (|104[) . as 
we shall explain.) The theory was worked out in full generality in [TS], [IH], [10] ■ 

Specifically, one starts with the Calderon formula: if c g (0, oo) is defined by c = Jq°° \ f{i)\^Y^ then for 
all u > 0, 

r°° df 

I \f{tu)\'j^c<^. (4) 

Discretizing (|4]), if a > 1 is sufficiently close to 1, one obtains a special form of Daubechies' condition: for 
all u > 0, 

oo 

0<Aa< J2 P <5a <00, (5) 

j = -oo 
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where 



A, 



c 



(l-0(|(a-lf(log|a-l|)|), 



(6) 



a 



2|loga| 



c 



(l + 0(|(a-l)2(log|a-l|)|)). 



(7) 



a 



2|loga| 



(([6|) and ([7]) were proved in [16], Lemma 7.6. In particular, Ba/Aa converges nearly quadratically to 1 as 
a ^ 1. For example, Daubechies calculated that if f{u) = ue~" and a = 2^^^, then Ba/Aa = 1.0000 to 
four significant digits. 

Let P be the projection in L^{S^) onto the space of constant functions (the null space of A.) Suppose 
that a > 1 is sufficiently close to 1 that Aa and Ba are close (or that a, f have been chosen in such 
a manner that Ba = Aa). Then for an appropriate discrete set {xj^k}{j.k)Gi.x[i.Nj] on M, and certain 
weights Hj i^, the collection of {(pj^k ■= l^-j^k^ai ,Xj ^1 constitutes a wavelet frame for (/ — P)L'^[S'^). By 
this we mean that there exist "frame bounds' Q < A < B < oo, such that for any F E {I — P)L'^{S'^) the 
following holds: 



If the points {xj k} are selected properly, B can be made arbitrarily close to Ba and A can be made 
arbitrarily close to Aa- Because of the proximity of Aa to Ba, the frame is therefore a "nearly tight 
frame" . In passing from one scale to another, we use "dilations" in the sense that we are looking at the 
kernel of /(t^A) for different t, which corresponds to dilating the metric on the manifold. On the sphere, 
the wavelets at each scale j are all rotates of each other, up to constant multiples. (This is analogous to 
the usual situation in which all the wavelets at each scale are translates of each other.) 
Again, wt,x is given by ([3]). On the sphere 5*^, we have 



The most important cases to consider on the sphere are the case in which / has compact support away 
from (the "needlet-type" case, of Narcowich, Petrushev and Ward), and the case in which f(u) — u^e~^ 
for some integer r > 1 (the "Mexican needlet" case, as considered in [18], [H]). (Actually, in their work, 
Narcowich, Petrushev and Ward used P in place of l{l + 1) in (I105|) . but this is a minor distinction.) 
Needlet-type wavelets and Mexican needlets each have their own advantages. Needlet-type wavelets have 
these advantages: for appropriate /, a, f. and Xj,k, needlet-type wavelets are a tight frame on the 
sphere (i.e., A = B in ([5])), and the frame elements at non-adjacent scales are orthogonal. In fact, for 
appropriate needlet-type /, a, fij and Xj,k, one can arrange that Aa = Ba = A ~ B = I in the discussion 
above. In that case, the wavelets are called needlets. 

Mexican needlets have their own advantages. In [18j . |19j . an approximate formula is written down for 
them which can be used directly on the sphere. (This formula, which arises from computation of a 
Maclaurin series, has been checked numerically. It is work in progress, expected to be completed soon, 
to estimate the remainder terms in this Maclaurin series. In [18], [19' the formula is written down only 
for r = 1, but it can be readily generalized to general r. The formula indicates that Mexican needlets are 
quite analogous to the Mexican hat wavelet which is commonly used on the real line.) Mexican needlets 
have Gaussian decay at each scale, and they do not oscillate (for small r) . Thus they can be implemented 
directly on the sphere, which is desirable if there is missing data (such as the "sky cut" of the CMB). 
It would be worthwile to utilize both needlets and Mexican needlets in the analysis of CMB, and the 
results should be compared. In this article we focus on needlet-type situations. 

A key property of this kind of wavelet is its localization. Working in the general situation (general M, 
general /), one has: 




(8) 



OO 




(9) 



i=0 m=~l 
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for every pair of C°° differential operators X (in x) and Y (in y) on M, and for every nonnegative integer 
N, there exists c := cn,x,y sucli that for all t > and x,y ^M. 

l^^^*(^'^)l^^(draAF' 

where / := degX, J := degl", and d{x,y) denotes the geodesic distance from x to y. In the case where 
M is the sphere, / has compact support away from 0, and if one replaces A; — l{l + 1) by P in the 
formula for Kt, this was shown by Narcowich, Petrushev and Ward, in [37] and [38]. The general result 
was shown in [18j . 

The proof of PH]) in [TH] would in fact go through without change if, instead of Kt being the kernel of 
/(i^A), where A is the Laplace-Beltrami operator, one assumed that Kt was the kernel of f{t^D), where 
D is any smooth positive second-order elliptic differential operator, acting on C°°(M.). It is natural to 
conjecture then that something analogous holds if D is instead a smooth positive second-order elliptic 
differential operator, between sections of line bundles on M. Since A^ may be interpreted as being such 
an operator, it is natural to want to prove the following result. 

Theorem \6A] 

Let Kt be the kernel of f{t^As). Then: 

For every R,R' G 50(3), every pair of compact sets C Ur and J-ri C Uri, and every pair of C°° 
differential operators X (in x) on Uu' and Y (in y) on Ur, and for every nonnegative integer N , there 
exists c such that for all t > Q, all x & Tri and all y G !Fr, we have 

|™,,,,,(.,,)|<c^g^, (11) 

where I := degX and J :— degY . 

In this article, we only give the proof of Theorem 16.11 in the easier "needlet-type" case, where / has 
compact support away from the origin. The general case will be established elsewhere. 

Let us now explain what we mean by spin wavelets. For R G 5*0(3), and x G Ur, let us define 

WtxR = X! X! /(^^'^''*) sYlmRix) sYlra- (12) 
i>|s| m 

Then wtxR S Of (5^). Also note that wtxR.,R'{y) = Kt.R,w{x, y) (if x e Ur and y G Ur>); this generalizes 
the case s = of ([3]). Moreover, if F G L%S^), then for x £ Ur 

Pt,F,x,R ■■= {F, m^R) = if{t^As)F)R{x) - {(3t,p)R{x) (13) 

where we have set Pt p — {f{t^As)F). We call PtFx,R ^ ■^P*"- wavelet coefficient of F. If F = 
T.i>\s\ E™ «i™ ^ 'lI{S'^), then 

Pt,F,x,R = f{t^hs)ai^ sYi^r{x). (14) 

l>\s\ rri 

Let us now explain how, in analogy to the case s = 0, one can obtain a nearly tight frame from spin 
wavelets. Let P = P\s\,s be the projection onto the H|s|,s, the null space of A;^ (in Of). Suppose that 
a > 1 is sufHciently close to 1 that Aa and Ba (of ([6]) and ([7])) are close (or that a, f have been chosen in 
such a manner that Ba = Aa). We then claim that for an appropriate discrete set {xj.k}(j.k)&i.x[i.N ■] on 
5^, certain Rj^k with Xj.k G Ur-^, and certain weights /i^ j., the collection of {4>j^k '■— l^j,k'^ai ,Xj ^-Rj k\ 
constitutes a nearly tight wavelet frame for (/ — P)Ll{S'^). (Note that this property is independent 
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of the choice of Rj^k-) In that case, we will have that for some C, if -F S (/ — P)L^, then, in L^, 
F ~ C'^^ f,{F, (pj I. (note that this sum is independent of the choice of Rj,k)- The precise statement 
is Theorem 16.31 below, where we explain how the Xj^k and fij are to be chosen. Theorem 16.31 will be 
proved elsewhere i2T] : in this article, it is only used in section 8. 

One can certainly choose a, and / with compact support away from the origin, so that Aa = Ba = 1 
(recall ((S])). Such / were used in [37] and [S^- It is then natural to conjecture, that as in the spin 
case of [33 and [3^, one can choose the points Xj^k so that the {^J'j^k'^j^k} ^n exactly tight frame 
for (/ — P)Ll{S'^), for suitable weights j,. One should note, however, that even if such explicit points 
and weights could be found, it might not be practical for physicists to use them in CMB analysis, 
where measurements are invariably taken at the HEALPix points [24j . This, then, is a situation where it 
is very useful (and perhaps essential) to have the flexibility in choice of the Xj^k which Theorem l6.3l allows. 

In section 7, we begin to examine how our theory can be applied to cosmology. We look at random spin 
s fields G on 5'^. We say that G is isotropic (in law) if for every xi, . . . , xn € C//, the joint probability 
distribution of Gf{xi), . . . ,Gf{xN) is independent of i? e 50(3). Similarly, if G, iJ arc both random 
spin s fields, we say that G, H are jointly isotropic (in law) if for every xi, . . . , xat, yi, . . . , um G C//, the 
joint probabihty distribution of Gf (xi), . . . , Gf (xat), Gf (yi), . . . , Gf{yM) is independent of i? S 5*0(3). 
Say G, H are jointly isotropic, and denote their spin s spherical harmonic coefficients by aim = {G^, syim), 
hm = {H^, sYim)', these are random variables. In Theorem 17.21 we show that E{aimbi'm') — unless 
I — I' and m = m' (see [6] for the scalar case). Moreover E{airnbim) does not depend on m; we will 
denote it by Ci^g,h- 

In Theorem 17.31 we generalize a key result in [4] , to show that the spin wavelet coefficients of G, H are 
asymptotically uncorrelated, under mild hypotheses. That is, for x € Ur, we define the spin wavelet co- 
efficients PG,t,x,R =< G,wtxR. >, PH,t,x,B^ =< H,wtxR >, and we show that |Cor(/3c.,t,a:,fl> /3ff,t,y,i?,)l ^ 
as t — > 0"*". (Again we assume that we are in the "needlet-type" situation, where / has compact support 
away from 0.) 

Finally, in section 8, we offer a brief glimpse of the consequences of the considerations of section 7 for 
stochastic limits and cosmology. 

Say e > 0. We use a nearly tight frame, as provided by Theorem 16.31 Specifically, in that theorem, we 
fix a real / supported in [1/a^, a^] for which the Daubechies sum J2JL-oo Pi^?'''^) = 1 for all it > 0, so 
that Aa = Ba = 1. We then produce the {(pj j,} = {0^ j, with frame bounds B — 1 + e and A ^ 1 — e. 
Say G = ^ Aim sYim is an isotropic random spin s field, and let Ci = Ci^g,G (notation as in Theorem 
17. 3p . For j e Z, we also let 7^- = ^a-^<a^^ \i,<a^ Ci{2l + 1); this is a small quantity since 
J2j Ij G;(2/ + 1) = 3VarG. (In all sums here, I satisfies the implicit restriction I > \s\.) We let 

Pjk — Pjkt — We focus on the quadratic statistics 

k k l,m 

When considering Fj, we can even let e depend on j here. In Proposition [8?T1 we show that 

E{\f,-r,\)<ej^ (16) 

for all j . 

Since we can let e << 1 depend on j, by Tchebychev's inequality, the probability distribution of Tj is a 
small perturbation of that of Tj . 

We conclude, then, by briefly discussing the latter. We assume now that G is Gaussian, by which we 
mean that {^Gi{x) : x £ Uj} U {3G/(x) : x G Uj} is a Gaussian family. We shall also assume that 
G = ^ Aim sYim is involutive, by which we mean that Aim = ^;,-m for all m (so that, in particular, Aiq 
is real). This assumption is very natural to make in cosmology. Indeed, as we shall explain in section 4, 
Newman and Penrose explained how to write any random spin s field F as F = F^, + ii^M, where the 
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"electric" and "magnetic" parts and Fm are involutive. In cosmology one assumes that these fields 
are jointly isotropic. 

On the other hand, the assumption of Gaussianity is common to a large part of the literature on random 
fields, and it is predominant in the cosmological area, both within physical and mathematical articles. 
Concerning the physical literature, Gaussianity is motivated by the widely dominating inflationary model 
for the Big Bang dynamics, which predicts the polarization random field is the outcome of Gaussian 
fluctuations of a quantum mechanical origin. From the mathematical point of view, Gaussianity is nearly 
a mandatory assumption in the high-frequency asymptotics framework we shall entertain here (compare 
[32]): indeed, as a consequence of the characterization of isotropic random fields given in ([6 ), the 
Gaussianity of the spherical harmonic coefficients is equivalent to their stochastic independence, which in 
turns make limit theorems feasible. The extension to non-Gaussian circumstances is certainly a crucial 
issue for further research; we stress, however, that in the framework of high-resolution asymptotics no 
rigorous results are known in the literature in non-Gaussian circumstances, even in the standard case of 
scalar random fields. 

Under these assumptions, in Proposition 18.21 we prove a central limit theorem for the Tj, namely that 




d N{0, 1) as j -oo , 



-^d denoting as usual convergence in probability law, and iV(0, 1) denoting as usual the standard normal 
distribution. 

Let us explain the significance of this result. As we explained earlier, our main motivating rationale is 
provided by the observations of a polarization random field F in cosmological satellite experiments. We 
shall then consider the realistic situation where a physical model is adopted to provide an expected value 
for the "angular power spectrum" C;, of the random spin field F, typically as a function C/ = C;(^'), of 
some vector of physical parameters ^. In the case of cosmological polarization data, the function C;(.) 
is estimated by the numerical solution of partial differential (thermodynamic) equations representing 
the evolution of a matter/energy fluid, and the parameters ^' are related to constants of fundamental 
physics. Propositions 18.11 and 18.21 highlight the possibility of using spin quantities, such as the spin 
wavelet coefficients /S^j, for F, to provide asymptotically valid estimators for scalar quantities of physical 
interest such as C/. For instance, by Propositions 18.11 and 18.21 it is immediate to construct a test for the 
assumption that the angular power spectra take the form C/(5'), for some physically motivated set of 
parameters which takes a given value 5* = ^!q. It is indeed sufficient to focus on 

_ Tj - F^qFj 

where E^^Tj, Var^g l^-^} expected mean and variance evaluated under the assumption that the 

vector of parameters is set equal to vPq. Propositions 18.11 and 18.21 predict that, for —j sufficiently large 
and e sufficiently small (depending on j), the probability of exceeding any given threshold value will be 
nearly provided by 

PrjS'j > Za/2} ^ a , 

where < a < 1 and Za/2 is the well-known quantile of the standard Gaussian distribution $(.), i.e. 
Za/2 '■— *I'~^(1 — ck/2)- Values of Sj above Za/2 provide evidence at significance level a against the as- 
sumption that represent the physical parameters in the model that generated the observations. 

We stress that the direct use of Tj is only feasible when observations from the whole sky are available. In 
cosmology, one is hampered by the presence of unobserved sky regions, for instance due to the foreground 
emission by the Milky Way and other astrophysical sources. The localization properties of needlet-type 
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spin wavelets in the real and harmonic domains are then clearly most valuable because, outside the 
masked regions, wavelet coefficients in the high-frequency limit are essentially unaffected by the missing 
observations. Moreover, the uncorrelation properties of needlet coefficients in the Gaussian case make a 
sound asymptotic statistical theory feasible, even in the presence of unobserved regions. Furthermore, 
these dual localization properties make possible the search for features and cosmological asymmetries, as 
done in the scalar case by [42]. We refer to [T^ for more discussion and details. (The reader should be 
aware that in the notation of 117], the high-frequency limit corresponds to letting j ^ oo, not — cx3 as in 
the present article.) 

3 Spin s line bundles 

In the next section, we shall review the groundbreaking work of Newman and Penrose |39| on spin 
weight s. They were writing for physicists, so it is necessary to place their observations in a rigorous 
mathematical context. We do this in a manner which remains very close to the spirit of the work of 
Newman and Penrose, and uses very little machinery; for another approach, using more machinery, see 
|13j . We also contribute some new results, some of which will be crucial for our work. 
Newman and Penrose work on a two-dimensional surface. Writing for physicists, they say that a quantity 
r] has spin weight s G Z, provided that, whenever a tangent vector to at a point transforms under 
coordinate change by to' = e^'^m, the quantity rj, at that point, transforms by 77' — e'^'^jy. Our first task 
is to put this notion into acceptable mathematical language. 

Let M be an oriented Riemann surface. The complex structure on M determines a unique metric (inner 
product on real tangent vectors at each point), up to conformal equivalence (that is, up to a constant 
multiple at each point). (Indeed, locally, such a metric is the puUback through a chart of the standard 
metric on C; different choices of charts lead to conformally equivalent metrics, since the Jacobian of a 
conformal mapping is always a constant multiple of an orthogonal transformation.) We always use a 
metric on M in this conformal class (that is, a metric which agrees with one determined by the complex 
structure, up to a constant multiple at each point) without further comment. Let us call this conformal 
class the intrinsic conformal class on M. 

li p G M, as usual we let Mp denote the real tangent space at p. Since M is oriented, for each v e Mp 
we may select Jv £ Mp at that point, such that {v, Jv) is an oriented orthogonal basis for Mp. Surely 
Jv is unique, up to a positive multiple. Thus, if v, w are nonzero tangent vectors at a point, and g is our 
metric, we may very naturally define the angle from v to w to be the angle -0 G ("''"i^r] determined by 
the two properties: 

(1) g{v,w) = rcos'0 for some r > 0; and 

(2) sgn V' = sgn g{Jv,w). 

Now, say we have an atlas U on M.. Suppose that for each chart Ua & U we have a section of 
TUa (the real tangent bundle of Ua). (This section will give, at each point of Ua, a "reference direc- 
tion".) Set p = {Pa}ua,eu- Say also s € Z. We then define L*, the spin s line bundle associated to 
(M^U, p), as follows. Say Ua, Up are charts in and define, for each p £ UaC] Up, ipppa to be the angle 
from Pa{p) to Pp{p)- We then set 

A/3o(p)=e*^^p^<^. 

Then XapXpa = 1 on n Up, and if U^ is a chart in U as well, then Xa-yX~/pXpa — 1 on Ua Ci Up U-y. 
By a standard argument (see, for instance, the argument on [26', pages 139-140), the Xpa may be used 
as transition functions to define a complex (not necessarily holomorphic) line bundle L'*; this bundle is 
unique up to isomorphism. If tt is the projection in L'* onto M, there are therefore difleomorphisms 
$Q : T^~'^{Ua) Ua 'X C, such that if p e L/q n Up, and if 
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then 

Let us say that these maps are induced by {U, p). 

Note that, if 17 is a open subset of M, we may naturally restrict L' to become a line bundle Lq over O, 
by using the charts {Uad^ U £U}, and the section p'^ ■= Pa J/a n f2. 

Suppose now that Lf, L| are determined, respectively, by (M\ZY\p^) and (M^, W^, p^). Say F : — > 

is holomorphic, and is a local diffeomorphism. Then F naturally gives rise to a smooth map Fg : lif 
L| as follows. Say (U , p^) induces the maps {$;!>,}, and {U , p^) induces the maps {^^j}- Say p e U^, a 
chart in , and q = F{p) G [/|, a chart in U^. The map F gives rise to a mapping F^p : M^. Let 

ip be the angle from i^*p[p^(p)] p^{q)- Say now {p,r]) G 7r~^(f/^), and that 

Then we define Fs{p,T]) = (g, C)j where 

Since F* preserves angles, it is easy to see that this definition is independent of choice of a, f3. Further, 
if / : Ml Ml is the identity map, so is Is : L^ L2. 

Moreover, say F is actually a biholomorphic map. Then, since and preserve angles, one easily 

sees that {Fs)~^ = {F~^)s. If / is a (continuous) section of Lj, we may naturally define a section of 
Lf by 

f^ = F-'ofoF (17) 
If G : M" M^ is another biholomorphic map, we then have 

irf = /^°^. (18) 



This is all quite abstract, but in fact we shall study spin s bundles only when the manifold is an open 
subset of the sphere 5^ or the complex plane C. On 5^, by choosing the atlas and p well, we can impose 

a very interesting and very explicit structure. 

On S"^, realized as {{x, y, z) € : x"^ + y"^ + z'^ = 1} , we let N be the north pole (0, 0, 1), and let S be 
the south pole, (0,0, —1). We define the chart Uj to be 5^ \ {N, S}. (The choice of chart map from Uj 
to C is irrelevant; one can take it to be a stereographic projection.) We obtain all other charts in our 
atlas by rotating Uj. Thus, if i? G 50(3), we define Ur = RUj. On Uj we often use standard spherical 
coordinates {0,cj)), and at each point p of Uj we let Pj{p) be the unit tangent vector at p which is tangent 
to the circle 9 = constant. (For definiteness, if this circle is a; = r cos cfi, y = rsuKp, z = constant, let us 
choose pj{p) to point in the direction in which </> increases, the "counterclockwise" direction.) On any 
other chart Ur, we choose 

Pr{Rp) = R.p[pi{p)]- (19) 

for any p ^ Uq. 

On the chart Ur we can use coordinates [Or, <I)r), obtained by rotation of the {9, (p) coordinates. Thus, if 
the {6, (j)) coordinates of p G Uj are (^0, <f>o), then the {9r, (I)r) coordinates of Rp G Ur are also (^0, <t>o)- 
Then, at any q G Ur, PR{q) is the unit tangent vector to the circle = constant, pointing in the 
direction of increasing 0^. Again, the angle ippR-^Rx is the angle from Pr^{p) to Pr^{p)- We would clearly 
obtain the same angle if Pr{p) had been chosen as the unit vector pointing in the direction of increasing 
9r, or in the direction of any fixed linear combination of the unit vectors pointing in the directions of 
increasing (pj^ and increasing 9r. In this precise sense, ^pR^R^ effectively measures the angle by which the 
tangent plane at p is rotated if one uses the {Or^ , (I>r2) coordinates instead of the {9r^ , <pR^ ) coordinates. 
This is why one should think of our choice of Pr{p) only as a convenient "reference direction". 
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On C we will use only a single chart, C itself, and if p g C, we let our reference direction p{p) be d/dy. 
Since there is only one chart, this bundle is trivial. However, if F maps an open subset J7 of C confor- 
mally onto another open subset U of C, we can consider the map Fs^ which, as we shall see, can be quite 
interesting. 

Now, if / is a (continuous) section of L'* over an open set f2 C S*^, then we may write 

for all q € Ur n fi, for a suitable function Jr : URni} ^ C Note, then, that if p e Ui H Ur n il, then 

fRip) = e^'^fiip) (20) 

where 

i(j is the angle from pj{p) to Pr{p) — R^.p{pj{R^^p)). (21) 

Note also that if C 5^, and Ra E 5*0(3), any smooth function f(Rg) : URg nil C determines a unique 
section / of L"* over URg n V, with /r^ = f{Ro). Indeed, say R £ 50(3). If p S Ur^ n n fi, then /rXp) 
must be given by 

fRip)^e^'^-fRgip). (22) 

where of course, /r^ = f(Ro)i and where 

tp is the angle from PRg{p) to Pr{p)- (23) 
Conversely, there is surely a section / of over Ur^ n Vl with these Jr. 
Consequently the space of smooth sections f of over may be identified with the space 

{/ = (/fl)i?GSO(3) : each e C^{URnn), and ^ hold for all Rq, R G 5*0(3) and all p G UR„nURnn}. 

(24) 

We shall frequently make this identification, which often greatly simplifies matters conceptually. However, 
as is often the case in mathematics, certain properties are clearer, and more easily verified, if one uses 
the coordinate-free line bundle point of view. 

Equation (|20p describes what happens when one "rotates the charts" , and we now argue that "you get 
the same answer if you rotate the section instead". Thus, if i? £ 50(3), as we have seen, R : ^ 
gives rise to a smooth map Rs : h'^ ^ . If / is a section of L'' over C 5^ , using ([17]) we obtain the 
"rotated section" 

/^ = i?;io/oi?. (25) 

over R-^n. (Note that if s = 0, /-^ = / o R.) Thus, for p e fl, 

(/«)(i?-V)=i?;'°/(p), 

from which we see that 

(/^),(i?-V) = e-V,(p) (26) 

where 

ip is the angle from R^p{pjip)) to pj{R^^p). (27) 
Note that this ip is the same as the ip in ([2T|) . since preserves angles. Therefore 

fRip) - (/^)/(i?-V). (28) 

One could take the point of view of ([24]) here, and simply define f^ by setting {f^)i{q) = fR{Rq). 
But then it would be tedious (though certainly possible) to verify that the smoothness of / implies the 
smoothness of or that (/^^)^^ ~ /^^^^, both of which are evident from ([25]) . 
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We also need to discuss the metric and the orientation we will use on S'^. Since rotations are holomorphic 
maps on S'^, any metric in the intrinsic conformal class on is rotationally invariant, up to a constant 
multiple at each point. Thus there is a truly rotationally invariant metric in this conformal class, ob- 
tained by rotation of the metric at any particular point on the sphere. We note that such a rotationally 
invariant metric coincides (up to a constant multiple) with the Euclidean metric, as restricted to tangent 
vectors on the sphere. (That is, for some positive constant c, if w = 'I2j=i^j'§^j ^^^^ ^ ~ '^j=i^j'§^ j 
are real tangent vectors at any point on the sphere, their inner product there is c < v,w >e, where 
< v,w >E= 2j=i Indeed, since < v,w >e is rotationally invariant, this fact need only be checked 

at the South pole S. We can determine a metric in the intrinsic conformal class on S'^, by pulling back 
the standard metric on C through stereographic projection from the North Pole onto the tangent plane 
to the sphere at S; denote this stereographic projection by T. It is geometrically evident, and it is easily 
calculated (see the formula ((5T|) below for a = T/2), that T*s(9/(?x) = d/dx, and T^s{d/dy) = d/dy. 
Thus the tangent vectors ^ and are, up to a choice of conformal factor at S, an orthonormal basis 
at S, as desired. We use < v,w >e as our metric on S'^. We choose the orientation on which makes 
T orientation-preserving. Thus {d/dx, d/dy) is an oriented orthonormal basis at S, and, consequently, 
{d/dy, d/dx) is an oriented orthonormal basis at N. 

In our first new result about the spin s bundle over S'^, we take a section / of L*, and examine the 
behavior of // as we approach the north and south poles. 

Theorem 3.1 Let f be a continuous section of L** over S^. Then the function e^'^'^ fi{9, (j)) extends 
continuously from Ui to Uj U and the function e^^'^'^ fi{9, (j)) extends continuously from Uj to Uj U S. 

Proof Let us show the second statement; the proof of the first statement is almost identical. For any 
(f> G [0, 2tt), let fi^(j,{9) — fi{6, ((>). It will suffice to show that limg^^- e~*'"^//^0(0) exists, uniformly in 0, 
and that this limit is independent of 0. 

For (f) e [0, 2tt), let 7^ be the great arc j^{9) — {9, (/>) {Q < 9 < tt). If p = -fcj,{9) is on this great arc, then 
Pi{p) points in the direction ^'^{9) x p (cross product); thus 

lim_p^(7^(6')) 

U — ^TT 

exists, and this limit is uniform in (j) (indeed, a rotation about N takes any one of these limiting situations 
into any other). Moreover, evidently, is a tangent vector at S, and given our choice of orientation at 
S, 

the angle from wo to is 0. (29) 
Now, let R be any rotation of which takes S to a point other than N or S. We then have that, for 

{9, ,p)GUin Ur, 

/7,0(^) = e-^^^«-*/fl,(0,0) 

where V'e ^ is the angle from Pi{p) to Pji{p), \i p — {9,4)). For fixed the limit &s 9 tt^ of the right 
side exists (uniformly in 0) and equals e~*'**'^/fl;(S), where now 

is the angle from to pji{S). (30) 

Thus for any (j), 

lim e-''*'^//.^(6l) = e-''*("*+*'^)/i?(S) =e-*"('^+*'^-*°)e-*'**°/i?,(S) =e-*'*('^+*"*-*«) lim fi.o{9) = lim // 

since VP^ — V^o equals — </> (mod 27r), by (f29l) and ([30]) . This proves the second statement; the proof of the 
first statement is entirely similar, if one takes into account the difference in orientation at N. 
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4 The Newman-Penrose Theory 

In this section we review the spin s theory of Newman and Penrose [39j . We do so within the rigorous 
framework we have presented in section 3. Most of our arguments are very close to theirs. We resume 
presenting new results in the next section. 

Although we earlier used stereographic projection T from the north pole onto the tangent plane at S, 
in what follows we will follow Newman and Penrose and instead use stereographic projection a from 
the north pole onto the equatorial plane, from [7/ C to C* := C \ {0}. Note that a = T/2 is still 
orientation-preserving. Explicitly, if p G C// has Euclidean coordinates {x,y,z), then 

crp= -T—ix + iy), (31) 
1 — z 

which one easily visualizes by using similar triangles. From this, and the formulas sin6' = 2 sin | cos |, 
1 — cos 6' = 2sin^ |, it follows that if p € J7/ has standard spherical coordinates {9, 0), then 

crp = e*'^cot(^). (32) 

The map tr, being a conformal map, preserves angles; we will also need to know the factor by which it 
multiplies lengths. Thus, if p G Ui, since a is conformal, there is a number Ap > such that for any 

1/2 

tangent vector v at p, \cr^pv\ = Xp < v,v >^ . We need to find the conformal factor Ap. For this, 
we may assume that v = Pi{p) is a unit tangent vector tangent to the circle <j) — constant through p. 
Restricted to that circle, by (j3ip . the map a is simply a dilation by the factor (1 — z)^^ , so we must have 
Ap = (1 — z)^^. It is easy to calculate, then, from ([3T|) . that Ap = 5(1 + lo'Pp)- For this reason, for C G C, 
we set 

^(C) = ^(1 + ICP), (33) 

and we have shown that: 

1 /2 

if p G C//, ap = C, and u is a tangent vector at p, then \a<,pv\ = P{() < v,v >^ . (34) 

We also remark the following fact (which is also geometrically evident). Since a restricted to the afore- 
mentioned circle is just a dilation, if 

li p = {6, 4>) G Ui, then the angle from (J*pj{p) to d/dy is — (j). (35) 

Indeed, it equals the angle from pj{p) to d/dy in M^, which is surely —(j>. 

We are now ready for the observations of Newman and Penrose, expressed in our language. These ob- 
servations will be preceded by bullet points in what follows. We will provide rigorous proofs of their 
observations, within the context of spin s line bundles. Our proofs are modified versions of their argu- 
ments. 

Say n C Ui is open. Suppose /i(7) : — > C is smooth. Following Newman and Penrose we define 

(Here we are again using standard spherical coordinates (0, 0) on [//.) Thus dsih^j^ is a smooth function 
on SI. 

Next, if is a smooth section of Lfj, we define dh to be the section of L''+^ over O such that 



{dh)i = dsihi. 
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Also, if Qq C C* is open, and m is a smooth section of Lqqj "^"^ define the smooth section u of L^^^ by 

(Here P is as in (133^ . Recall that we have only one chart on C*, so sections of Lq^ or L^^^ may be 
naturally identified with smooth functions from flo to C.) 

If / is a section of L^, where C Ui, let us write ) — f °s <J~^ ■ The first key observation of Newman 
and Penrose is: 



Suppose / is a smooth section of , where C [// is open. Then 



(djy -i?,(r" ). (36) 

To see this, one notes that dsiji — scotOfi — Lfj, where L = {d/d9) + [i/ sin6){d / dcf)). Using ([55l) . 
we see that the left side of ^ is e-^^'+'^'i'^[{s cot 9) fi o cj-i - (Lfi) o a^^]. Let m = // o ct"^ A 
brief calculation, using ([5^ . the chain rule, and the fact that cot{9/2)/ sm9 = icsc^(0/2), shows that 
(Lfi) o CT'i = -e">'csc^{9/2)du/dC. On the other hand, if p = (6',</>) e Uj and C = crp, then P(C) = 
i[l + cot2(6'/2)] = i csc2(6'/2). Thus the left side of 1^ equals 



-i(s+l) 



(scot6')u + 2Pe"^|^ 



On the other hand, the right side of is 



2P' 



dC ' 



which equals 2pi-"[5e-'""^PVaC]u + 2Pe-'">'du/ dC,. We are left with showing that 

^^_^de-^"f'P' 



9C 



scot( 



(37) 



(38) 



(39) 



Surely dP/dC, = \(^. On the other hand, since e^*-^ = C/C, we find 2ie^^'l'd(j)/dC = 1/C, so that = 
l/(2iC). Thus the left side of ^ equals 



e-'-'P[-sP/C + sC] = 



'[m-p)/c] = 



■[(|Cp-P)/cot(0/2)], 



so ah that is left to show is that (|CP - P)/cot{0/2) = cot9. But 
(ICP - P)/cot{e/2) = [cot^{9/2) - l]/[2cot(6l/2)] = cot 6* as desired. 

Next, Newman and Penrose observe: 

• Say R G 5*0(3). Suppose / is a smooth section of L^, where ft C Uj H Ur is open. Say R G 50(3). 
Then 

= {dff- (40) 
(Note that these sections are defined over R^^Vl QUi ) 

To see PO)) . it sufhces to show that [d{f^')Y — [{df)^Y ■ (Note that these sections are defined over 
rii := <j{R^^Vl).) Set u ~ , defined over Hq := ail, and let 



Then by (pS]). 



X = (7 o Ro a ^ : Qi ^ flo 



(41) 



15 



On the other hand, 



a oaoRoa 



{Vs uf. 



Thus, we need only show that if u is a section of Lf^^ C C*, then 



(42) 



(Note that these sections are defined over fii.) 

Now X is a conformal mapping. At any point C of fii, if !'{(,) ~ re'"'^ {ip G [— 7r,7r)), then wc know that 

the angle from l^(^{d/dy) to d/dy is — -0, (43) 
while for any real nonzero tangent vector v at — r\v\. By (j34p and (|4ip . we conclude that 



p(c) ■ 



Accordingly 



P(C)e 



_ Pirn) 



(44) 



(45) 



Turning now to ((42|) . we see by (l45|) that 

,i-.^[™e^'^)^(woT)(C)] _ 2P(C)^ 



(P, u^)(C)-2P(C)^ 



([P^u]oX)(C)/X'(C)' 



5C 



P(C)^ 



dC 



The key point is that T'{C,Y is antiholomorphic, so that 

2P(C) 9([P^zi]oX)(C) 2P(C) 



P(X(C))''e-^'"^ 



9(P"u) 



5C 



Finally, then, by (|44)) . 

(P, ti^)(C) - 2P(X(C))i-^e^(^+i)''' 
This proves (|42)) . as desired. 



(2?. ")^(C)- 



Next, say £7 C C/j^ is open, so that P ^£7 C [//. Say : ^ C is smooth; then /i(j^) o R : R ^17 — > C 
is smooth, so we can apply 9^/ to it and obtain another smooth fmiction on P^^i7. Now we define 

dsB^R) = [dsi{h(R) o R)] o R-\ (46) 
which is, like /i^^) itself, a smooth function on f7. We now deduce: 

• Suppose / is a smooth section of L" over an open set C C// n Ur. Then 

{df)R = dsRjR. (47) 

To see this, say p e 17. Then by ^ and (gOl), 

{df)R,{p) - [(a/)^],(p-V) = (a/^)/(p-V)- 

Thus, by the definition of dsi and ([28|) again, 

(5/)i?(p) = (a.z/f)(p-V) = (a.z[/flo i?])(i?-V) = {dsRjRm, 
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as desired. 

We can reformulate the last result more appealingly as follows. On the chart Ur we again use coor- 
dinates {Ofjjipj^), where if the {0,ip) coordinates of p g Uj are {9q,(I)q), then the {9ii,(j>ji) coordinates of 
Rp & Ur are also {9o, (/)q). Then, in these coordinates on Ur, 

dsnfn = -(sin..)^ (J- + -^^) (sin.,)-/.. (48) 

So far we have only defined d over the chart t//, but (|47l) enables us to define df for any smooth section 
/ of Lq, where now O is any open subset of . We simply let df be the smooth section of Lj^"^^, with 
the property that if p e is in Ur^ then 

{df)R{p)=dsRfR{p). (49) 

Note that this definition makes sense. Indeed, over n C//, this section agrees with the section df wc 
have previously defined. If, instead, p is the north or south pole, choose a particular R S 50(3) with 
p e Ur, and let h be the section of over il n Ur such that Hr — dsRfR- Then on n Ur DUi, h agrees 
with df as previously defined. Thus, for any other rotation R' such that p e Ur', hRi(q) = dsR' fR'{q) 
for all g e n Ur n Uri n Ui. Letting q ^ p, we see that hR'{p) = dsR'fR'{p) as well. This shows that 
the definition (j49|) makes sense. 
As a consequence we have that: 

d maps smooth sections of to smooth sections of L'*"'""'^. (50) 

Here the sections could be defined over the entire sphere; d is consequently called the "spin-raising 
operator" . 

Note also that, in (|40|) . we may now relax the hypothesis that C UiD Ur, and only assume that C S*^ 
is open. Indeed, both sides of agree on Q Ci Uj H Ur; since both sides are smooth sections of L^, 
they must agree on all of $7. 

As is customary for sections of line bundles, we let C°°(L^) denote the space of smooth sections of L*, 
defined over all of S'^. (In section 2, we called this space C^{S'^).) Newman and Penrose make the 
following important observation: 

• If s < 0, then d : C°°(L") ^ C°°{L'+^) is injective. (51) 

To see this, say / g C°°(L''), and that df ~ 0. By Theorem l3.1l |//| is bounded on Ui. Let u = // o cr^^; 
then u is bounded on C*. As in ([36|) and (l38|) . we see that 

= (df)"'' = V, (.r"') = 2P^-'[de-'"*'P'u/d(:]. 

The function v := e^'^'^P^u is therefore holomorphic on C*. It is bounded near 0, so it extends holomor- 
phically to C. But, since s < 0, [P(C)]* ^ as C — > oOj w = by Liouville's theorem. Consequently u, 
and hence /, are both zero, as desired. 

Say now that C S*^ is open. There is evidently a well-defined map of complex conjugation from 
C°°(Lq) to C°°(L^''), taking a section / in the former space to a section / in the latter space, defined 

by _ 

{f)Bip) = fnip) 

for all p e Ur n Vt (for any R e 50(3)). We may then define the "spin-lowering operator" d by 

df = 97, (52) 
so that d : 0°°(Lj2) 0°°(L^^^). In the coordinates {9r, on Ur, we evidently have that 

(df)R = dsRfR, 



17 



where 

-dsnfn = -(sm.,)- - (sin..r/.. (53) 

If Qo C C* is open, and it is a smooth section of L^^, we define the smooth section Vg u of L^"^ by 
Vs u — Vg u\ thus 

By taking complex conjugates in (I40p . and ([?T|) . we find: 

• Suppose / is a smooth section of Lq , where 17 C L// is open. Then 

[dfY'' =Vs{r-'). (54) 

• Suppose / is a smooth section of L^, where C 5^ is open. Say R <E 5*0(3). Then 

a(/«) = a/«. (55) 

• If s > 0, then a : C°°(L*) ^ C°°(L''-i) is injective. (56) 

One can define sections of L'*, (or on any other hne bundle over a compact oriented manifold, for that 
matter) , and define an norm on them, by using a finite atlas on the bundle and a partition of unity 
on the manifold. Different choices of atlas and partition of unity lead to the same space of sections 
with an equivalent L'^ norm. For the bundle L'* there is a very natural choice of norm, equivalent to 
the Li^ norm obtained through such a construction. Note first that if / is an section of L'', then (a 
representative of) / is defined almost everywhere. Say / is defined at p e Ur n Uri. We may then write 
^Ri.f{p)) = (pJrIp)), *fl'(/(p)) = {pJr'{p)), fR'ip) = e'^'^/i?(p), where V is the angle from p^(p) to 
Pri{p)- Thus |/fl,(p)P is independent of the choice of R, as long as p € Ur. Accordingly 

\f{p)?--^\fR{p)? (p^Ur) 
is well-defined a.e. on S*^. We then define the norm of / by 

\\f\\h= I \f? 



where the integral is with respect to the usual surface measure on S^. We let L'^{L^) denote the space of 

sections of L'*, with this norm. 
We may now show: 

• The formal adjoint of d : C°°(L'*) ^ C°°{L'+^) is -d. 

For this, we need to show that, whenever / e C°°(L*) and g E C°°(L''+^), we have 

<df,g>^~ < f,dg>, (57) 

where < •, • > denotes inner product. 

For this, we first note the following useful proposition: 

Proposition 4.1 (a) Let N be a positive integer. Then S'^ is not the union of N open subsets of diameter 
less than t:/N . 

(b) Any open set ft C_ of diameter less that tt/2 is contained in some Uu. 
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Proof For (a), say it were, and let Ui be one of those open subsets. By connectedness of S'^, Ui must 
intersect another one of those sets, say U2- Again by connectedness, Ui U U2 must intersect another one 
of those sets, say C/3. Continuing in this manner, we obtain the contradiction 

TT = diamS"^ < diamf/i + diamf/2 + . . . < Ntt/N. 

For (b), we need to show that omits some pair of antipodal points in S^. If it did not, then if A is the 
antipodal map, fl U Afl = S^. This contradicts (a). 

Let us then cover S*^ by a finite collection of open sets of diameter less than 7r/5, choose a partition 
of unity {(j} subordinate to this cover, and break up / = ^ Cj/ and g = J2 Cj9- By Proposition HH] (b), 
by using this breakup, in proving ()57|) . we may assume that supp/ U suppg Q Ur for some R. Then we 
need to show that 

27r pTT 



[{dsRfR)gR]{OR, (Pb) smeiid0Rd(f>jj = - 
But this is an elementary calculation, using (|48l) and ([SS]). 



[.fRds+i,RgR] {Or, sin OndORdcj)]^ 



We are now ready (again following Newman and Penrose) to define the spin s spherical harmonics. 
We let {Yira : < Z, — Z < TO < /} be the usual orthonormal basis of spherical harmonics for S*^, 



Yin 



■sn?\ei2)Y^{[) (-l)'-'-cot2-™(^/2). 



(58) 



where 

a;™ = [il + - m)\{2l + l)/(4^)]' , 
and our convention is that the binomial coefficient (") is zero if j < (or n < 0). (As usual to+ = 
max(TO, 0), TO^ — 'max{—m,Q). The factor of (—1)"' in (|58p is usually not included, but by including 
it, we have Yim = i1,-m, which will simplify many formulas in the sequel.) 
For I > \s\, Newman and Penrose then define the spin s spherical harmonic sYim by 





[(i-s)!] 


< 


[(i+sy.\ 




\{i+sy.] 


1. 


[{i-sy.\ 



d'Yirn (0 < S < I), 

{~d)-'Yhn {~1<S< 0). 



(59) 



The ^Yjm are not defined for I < \s\. Note that ^Yj„, is a smooth section of L'*. We have 



sYlml{9, I 



aw. 



^ ,m'\9/2) ^ (,|t-„0 (-l)'-"--^ cot^' 



2r-\-s — m 



{9/2). 



(60) 

To see this, let gVim denote the section of L'' over [// such that sVimi equals the right side of (pD|) . We 
need to show that sVim = sYim over Ui. Note that this is surely true when s = 0. Let gUim = svf^ ■ By 
(|36)) . when s > 0, we need only show that 



il + s)\ 



(Here V — I?s-i o ■ ■ ■ o 

sUlmiO = 



{i-sy. 

Vq.) However, using ([35]) again, one easily calculates that, whenever / > 



(61) 



-1)' 



at 



[{l + s)Kl-s)l 



J T- 



(62) 
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From this, a brief calculation using Pascal's rule and the identities 

(/ - . - r) i'r) = ir+l) a+i) = {l-s) 
shows that, whenever I > \s\, 

Vs sUlraiO = [{I -S){l + S+ 1)]^ s+lUlm- (63) 

(Here, our convention is, that since the factor (/ — s) appears on the right side of ([63|) . the right side is 
to be interpreted as zero if ^ = s.) From this, (|6ip follows at once by induction on s, for s > 0, and so we 
have (pO]) for s > as well. We then obtain (pOj) for s < by complex conjugation and use of the identity 



spiral —s^l. — 'tn.I- 

Even better, by and ([551) . '^^ have that, whenever I > \s\, 

dsYlrniC) = [{1 - + S + 1)]^ s+lYlrn- (64) 

By taking complex conjugates of this, we also find that, whenever I > |s|, 

dsYlmiC) -[il + s)il-S+l)]K-lYlrn. (65) 

Newman and Penrose next observe that, for each s, the {sYim} are orthonormal. Indeed, set 

_/[(/ + .)!/(/-.)!]V2 if <,,</, 



By (IMl), if s > 0, 



while if s < 0, 



If s > 0, we see that 



sY,m - ^d'Ylm, Y,m = 7^ ("5)^ sF/™, (67) 
Ols Ols 

sYi^ = ^{-d)-'Yi^, Yi^ = ^d-' ,F/„,. (68) 
his bis 



< sYlm, sYl'm' >— - — < d'^Yim, sYv m' >~ ; — < ^^m, { — dY sYl'm' > — < Yijn,Yif,n' >= Su'Smm', 
Ols Ols 

which proves the orthonormality. If s < 0, one obtains the orthonormality through complex conjugation. 
Thus, the {sYim} are an orthonormal set of eigenfuntions for the (formally) positive self-adjoint operator 
—dd : C°°(L*) C°°(L'*). Since this operator is elliptic (by ((^5)) and ([551) ). it has an orthonormal basis 
of smooth eigenfunctions, so at this point one would surely expect: 

• The { sYim} are an orthonormal basis of smooth eigenfunctions for — dd : C°°(L*) C°°(L'*). (69) 

This would follow if we knew: 

• If s > 0, then -.0°^ ^ C°°(L") is surjective. (70) 

Indeed, say s > 0. Now : C°° C°°(L*) is continuous (using their Frechet space topologies). (As 
usual, this topology is obtained by using a finite atlas on the bundle and a partition of unity on the 
manifold. Different choices of atlas and partition of unity lead to an equivalent Frechet space topology.) 
Thus, since the linear span of the {Yim} is dense in C°°, the linear span of the { sYim} is dense in 
Qs(jca _ (j°°{Jj^)^ which is in turn L^-dense in the space of sections. This would prove ([S^ for s > 
(and hence for all s by complex conjugation), if we knew ([70)1 . Moreover, ([70)1 itself follows from the 
injectivity of 9 : C°°(L") ^ C°° (see (HH)), and by use of elliptic theory and Sobolev spaces. 
Rather than proceed in this manner, and in order to extract important additional information, we give 
a direct proof of (|M)) and (|70p. by following the method of Newman and Penrose. Their method was 
somewhat formal, but it can easily be made rigorous. In fact, we shall show: 
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Theorem 4.2 Fix an integer s. 

(a) Let {II ||jv} be a family of norms defining the Frechet topology on C°°(L'*). Then for any N there 
exist C, M such that 

\\sYi„^\\n <Cil + l)^' (71) 

whenever I > \s\. 

(b) For every f G C°°(L*), there exist numbers {Aim ■ I > |s|, —l<m<l} such that 

; > I s I m 

(convergence in C°°(L'')j, such that the Aim decay rapidly, by which we mean that for every N there exist 
C, M such that 

|^/™| <C(/ + 1)-^^ (73) 
for all l,m. Further, we must have Aim =< /, sYim >■ 

Proof (a) is well-known when s = 0. It then follows for all s by use of ([59|) . as well as (|53|) and (|48|) for 
general R. Note, then, that any series as on the right side of ([7^ . where the Aim decay rapidly, does 
converge in C°°(L*)). 

(b) is also well-known when s = 0. To prove it in general, we may as usual assume s > 0. In that case 
we may write 

^'f^Y^Y. BlmYlm BlmYlm + 9 

i>0 rn l>s m 

where the B/m decay rapidly, and where g = J2i<s^m ^ kerc?''. Now, with bis as in set 

Aim = (—^y Bim/bis (for I > \s\). Then, again using (a) for s = 0, wc find 

^'(/ - EE -^'™) =9^ kcra^ (74) 

l>s m 

Now it is a simple general fact, that if D is a differential operator mapping smooth sections of one 
hermitian line bundle, over a compact oriented manifold, to smooth sections of another hermitian line 
bundle over that manifold, and if DF is in the kernel of D* for some smooth F, then DF = 0. (Indeed, 
D*DF = 0, so the inner product of DF with itself is zero.) Thus the left side of (|74p is zero, and so we 
have ([72|) . by ([56|) . Surely Aim =< /, sYim >, by the orthonormality of the sYim- This completes the 
proof. Moreover, ([70| follows as well, since if s > 0, 

l>s m 

if Clm — Aim /bis- 

As we have said, the completeness of the {sY/m} as an orthonormal basis follows from ([70)) . It also follows 
directly from Theorem l4.2l (b), which shows that the linear span of the { sYim} is dense in C°°(L*). 
If s > 0, we may write 



where 
and 

Moreover, 



V = rC°°(L^) = {Yim ■■l>s) 

W = kerd' ^ (Yim-. I <s). 
Qs .v-^C^iL") bijectively. 



(75) 
(76) 
(77) 
(78) 
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The statements ([75)1 and (|78l) are evident, once one identifies V and as spaces spanned by certain 
spherical harmonics, as in (|75)) and ([77]) . One could obtain ([75)1 and ([75]) by using only the first definitions 
of F and W in ([76]) and ((77)) . by using elliptic theory and ([SG]) . However, by identifying them as spaces 
spanned by certain spherical harmonics, Newman and Penrose observe a crucial additional structure. 
Namely, the decomposition (I75|) respects the splitting of a function into its real and imaginary parts. That 
is, ii g E V (or W), so are dig and ^g. (This is evident, since Yim — ¥/.-,„.) This splitting of any function 
in V then induces, through the isomorphism in (|78p . a corresponding splitting of any element of C°°(L*). 
Thus, say / G C°°(L*). Suppose s > 0. We write / = d'^g for some g G C°°, then set 

fE = d'ng, /m = 5^%. 

/e and /m are well-defined, independent of choice of g. Indeed, if d'^h = 0, then h G W, so 3?/i and 3/i 
are in W, so d'^dih = d'^^h — 0. In particular, we could let g be the unique solution of / = d'^g in V. 
Similarly, if s < 0, we write f — d g for some g S C°°, then set 

/e = d^^g, fm = d^^g. 

Surely we always have f — f^ + */m- 

Again suppose s > 0. By repeated use of ((55|) . one sees that the space V is rotationally invariant, and 
then by repeated use of (|40|) . one sees that 

{f"h^{fEf, (nM-(/M)^ (79) 

for any R E 5*0(3). Similarly these relations hold for s < as well. 

In fact, if / is as in (|72p. it is very easy to find the expansions of /e and /m in terms of the gYim- Say 
first s > 0. The solution g oi f = d'^g, with g E V, is simply 

l>s m 

A very brief calculation, using Y/^i = Y/.-m, now shows that 

/e- 5m^""E.lim (80) 

/>|s| m 

and 

/m = ^ ^ ^/?riM s^m (81) 

/>|s| m 

where 

A;™e = (^/m + v4;,_™)/2, (82) 

while 

^ZmM = -i(^/m - Ai,_™)/2. (83) 

Similarly, ([50]) . ([5T|) . ([5^ and ([55]) hold for s < as well. Up to normalization constants, these formulas 
agree with formulas in section II of 53J, and also with equation (54) of _7J (to see the latter, use equations 
(55), (65), (66) and (69) of that article). (Those articles take s — 2 and write f = Q + iU. |7j uses the 
notation a^^^ and a^^^^^ in place of our A;„je and A;,„m-) 

Let us now let 

' ^ I -dd if s < 0. ^ ' 
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(Here Ag acts on smooth sections of L'*.) We also let 



We also let 
Also 



^ ^Ul-s){l + s + l) ifO<s<Z, 

|(Z + s)(;-s + l) if - l<s<0. ^ ' 



His =<sYim:-l<m<l> (for I > \s\). (86) 



Pis denotes the projection in L^(L''') onto His- (87) 

Then L'^(L'^) — ©;>|s| Ti-is- We also have that His is the space of smooth sections of L* which are eigen- 
functions of for the eigenvalue Xis- (This follows from (f64|) . ([65]) . and Theorem 14.21 ) (In particular, 
then, Ao is the usual spherical Laplacian on S'^.) This was all noted by Newman and Penrose. 

Of course {sY/m : ~l < m < 1} is an orthonormal basis for His- Note that the map f ^ ^ — RsofoR^^ 
is a unitary mapping of L^(L*) onto itself, for any R € 50(3). Note also that for each such R, sYim 
is an eigenfunction of A^ for the eigenvalue Xis- (This follows from ([40|) and ([55]) .) Thus, for any R, 
{s^im ■ —I < m < 1} is also an orthonormal basis of His- Note also that, by 1^5]). we have, for p e Uji, 

{sYi'^")r{p) = sYi^i{R-^p)- (88) 
and so, in the {Br, (f>jf) coordinates on Ur, we have 

isYiV)RiOR,^n) = e''-'^H-ir' n, ,uT M11/2 E W Ut-™) (-l)'-"-^ cot2'-+^-™( 



[(; + s)!(;-s)!]i/2 



(89) 



5 Zonal Harmonics and Kernels 

We now return to our new results. 

Concerning the sYim, we also make the following observation, which is essentially in [7], section 6. 
Proposition 5.1 

s^a+ 10 otherwise 

Proof Since sYimi{S,4>) equals e*™"^ times a function of 9, this limit (which surely exists, by Theorem 
13. ip must be zero if to ^ — s. If to = —s, it follows from an examination of (j58p . (In each term in 
the summation in ([60]), because of the binomial coefficients, we must have r < I — s- As 6* ^ O"*", 
sin2'(6l/2) cot2'^+''-'"(6l/2) = sin2'(6i/2) cot2''+2''(6'/2) has a nonzero limit only if 21 < 2r + 2s. These two 
conditions together allow us to look only at the term in which r = I — s when taking the limit, and the 
proposition now follows at once.) 

We next discuss the kernel of the projection operator Pis] we claim that it is a smooth-kernel operator. 
Since Pis maps sections in L^(L*) to smooth sections of L*, we must define what we mean by the kernel 
of such an operator. For x € S"^, let be 7t~^{x), the fiber above x. 

Definition 5.2 We say that T : L^(L'') C°°(L^) is a smooth-kernel operator ifT has the form 



{Tf)(x)= f K(x,y)f{y)dS{y) 



(90) 
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for any f G L^(L'*); here we require: 

(i) For each {x,y) € x 5^, K{x,y) : — > is linear; 

(ii) K{x,y) depends smoothly on {x,y) G S"^ x 5^. 

We also call the restriction of T to C°°(L*) a smooth-kernel operator. We call the collection of linear 
maps K = {K{x,y) : {x,y) G 5"^ x S"^} the kernel ofT. 

Note here that the integral in ((90| is a vector-valued integral. Also note that the statement {ii) means 
that the kernel is smooth in (x, y) after one locally trivializes the bundle. 

To say that T is a smooth-kernel operator is evidently equivalent to saying the following: for each 
i?, R' G 5*0(3), there exists a function Kri^r G C°°{Ur' x Ur), such that for all smooth sections / of 
with compact support in Ur, we have 

{Tf)R,{x)^ ( KR,^R{x,y)fR{y)dS{y), (91) 

for all X G Uri. 

Evidently, if i?, i?', Ri,R'^ G 5*0(3), if ^ holds, we must have 

KR,>^R,{x,y) = e^^'^^-^^^KR^Mx.y) (92) 
for all X G Ur^i n Uri and y G Ur^ n Ur, where 

tp^ is the angle from Pri{x) to Pr-^i{x), and -i/jj is the angle from PR{y) to PR-^{y). (93) 
Let us now examine the kernel of the projection operator Pig. Evidently, if / G L^(L*), we have 

I 

Plsf = ^ < /, sYlrr. > sYl„^. (94) 

m— — l 

R-l 

(Here one could use si^„° in place of gYim, for any Rq G 50(3).) 
Thus Pis has kernel K'^", where for any R, R' G 50(3), 

I 

KR'^R{x,y)^ sYiraR'{x)sYi^R{y). (95) 

m—~l 

(Again here, one could use s^jm" in place of sYim, for any Rq G 50(3).) 
Definition 5.3 For \l\ > s, we define sZi, the "s-zonal harmonic" in Ti-is, by 

sZi = {-ly 

We have the following new theorem, which provides justification for this definition. When s = this 
theorem reduces to the well-known result that 

oZi{x) := Zi{x) = K^^{x, N) is the zonal harmonic in TiiQ. (97) 



2/ + 1 

An 



.Y, 



(96) 
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Theorem 5.4 (a) Say N e Ur, and that x e Uj n Uri. Then 

,Zu{x) = e'^('^i-'^^)ifj|,,^(x,N), (98) 

where 

ipi is the angle from Pjii{x) to Pi{x), and -02 *s the angle from Pn^N) to d/dy (at N). (99) 
(b) For any f G His, we have 

lim e'^'^f 1(9, cj,) =< f, sZi> . (100) 

9 — >0 

Proof For (a), note that if y e Uj Ci Ur, then K'j''j{x, y) equals 

J2 sYi„ri{x)sYi,r.i{y) = e''^^^-^-^^K'^,^R{x, y) 

m 

where ipi is the angle from Pr,{x) to pj{x), and t/ij is the angle from p^ijj) to Pi{y)- In this equation, 
let us say y has coordinates (6',(/)) (in the usual coordinates on [//). Let us fix 0, multiply both sides 
of the equation by e"'**^, and take the limit as 6* — > 0+. (One can, if one wishes, choose to take = 0; 
or one can let (j) be arbitrary, and recall our choice of orientation at N.) (a) now follows at once from 
Proposition 15.11 

To prove (b), one may assume / = sYim for some m; but then the result is evident from Proposition 
15.11 (Alternatively, one can start with the equation fi{x) = J^a K'fj{x,y)fi{y)dS{y) (which is valid for 
/ e His, X E Ui), multiply by e**"^, and take a limit.) This completes the proof. 

For another new result, also familiar when s = 0, we examine i^j|^(x,x), for x G Ur. It is evident 
from ((92|) . ([93]) that this quantity is independent of R. In fact it is independent of x as well: 



Theorem 5.5 If x eUr, then 

91 -I- 1 

K^AM^,x) = ^. (101) 

Proof It is enough to prove this for R — I, since we will then know it unless a; = N or S, and the result 
will follow in those cases too by continuity. Let us then take R = I and free the letter R for other uses. 
Say now x,y E Uj. Choose R G 50(3) with Ry = x. Then x £ Uj H Ur, so 

KYj{x,x) = K%r{x,x) = J2 \isYil'')Rix)\' = J2 Uyimi{R-'x)\' = K\^j{y,y). 

m m 

(We used Accordingly K\^j{x, x) is independent of a: € 50(3). Thus for any fixed x G 5*0(3), 

AtiKY,{x,x)^ [ K'/j{y,y)dS{y)^y2 [ \sYUy)\^dS{y) = 21 + 1, 

as desired. 

In doing estimates involving the if'**, it will be very useful to observe the theorem which follows. In 
this theorem, d^^ will denote the differential operator on Ur which satisfies {d'^F)R — O^^Fr, for any 
ordinary smooth function F on Ur; explicitly, d^^ — ds-\,R ° ■ ■ -da^R. Similarly, d^R will denote the 

—s —[s] 

diff'erential operator on Ur, which satisfies {d F)r — Or Fr, for any ordinary smooth function F on Ur. 
Theorem 5.6 Say x G Ur' and y G Ur. Then if s > 0, 

K'A'.ui^^v) = Jf^/wjR,yK'\^.vl (102) 
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while if s < 0, 

K%A^,y) = (103) 

(Here 9^] ^ denotes 9^] taken in the x variable, etc.) 
Proof This follows at once from (leTj) . ((68l) . and ([95l) . 

Of course K^^{x^y) = Zi^x{y), the zonal harmonic in TLiq based at x. It equals [{21 + 1) / A'K]P^^'^ {x ■ y), 
where P;^ is the ultraspherical (or Gegenbauer) polynomial of degree I associated with 1/2. 
Series defining wavelets on the sphere, involving the i^T'^, were estimated in [37], [3H] and [H]. Using 
these estimates, and Theorem 15.61 we will be able to estimate similar series, which involve the i^T'* for 
general s, in place of the iiT'", in the next section. 



6 Spin Wavelets 

In sections 3-5 we have often used the variable "/" to denote a spin function on the sphere. In the 
remainder of this article, the variable "/" will be reserved for another purpose. 

Specifically, say / € 5(]R+), / 7^ 0, and /(O) = 0. Let (M,(7) be a smooth compact oriented Riemannian 
manifold, and let A be the Laplace-Beltrami operator on M (for instance, the spherical Laplacian if M 
is 5^). Let Kt be the kernel of /(i^A). Then, as we have explained in section 2, the functions 

WtAy)^Ktix,y), (104) 

if multiplied by appropriate weights, can be used as wavelets on M. In case M = S"" and / has compact 
support away from the origin, we shall say these wavelets are "needlet-type" . 
On the sphere 5^, we have 

00 / 

Kt{x,y)=Y, E fit'^l)ylMYl^^{y)■ (105) 

1=0 m=-l 

As we have explained in section 2, a key property of this kind of wavelet is its localization. Working in 
the general situation (general M, general /), one has: 

for every pair of C°° differential operators X (in x) and Y (in y) on M, and for every nonnegative integer 
N, there exists c := cn,x,y such that for alH > and x,y ^M. 

f-{n+I+J) 

l™*(^'^)l^^(^raAF' ^'''^ 

where / :— degX and J :— degF. In the case where M is the sphere, / has compact support away from 
0, and if one replaces A; ~ l{l + 1) by P in the formula for Kt, this was shown by Narcowich, Petrushev 
and Ward, in [37j and [S^- The general result was shown in [18]. It will be important to note that in 
()106p . the constant c can be chosen to be < C||/||m, for some M. (Here we have chosen a nondecreasing 
family of norms {|| to define the Frechet space topology on 5o(K+) = {/ e 5(M+) : /(O) = 0}.) 
Indeed, this fact may be read off from the proof of (|106p . More simply, one may use the closed graph 
theorem: The collection of functions Kt{x,y) (smooth in t,x,y for t > 0) which satisfy (|106[) is itself 
naturally a Frechet space, say F. We want to see that the map taking / to the function Ktj{x, y) E F is 
continuous. (Here we are temporarily letting Ktj be the kernel of /(t^A).) For this, by the closed graph 
theorem, one need only note that the map from 5o(R^) to C which takes / to Kt,f{x, y) (for any fixed 
t,x,y) is continuous; but this is evident from ([TUSj). 

The proof of (|106p in [18] would in fact go through without change if, instead of Kt being the kernel of 
/(t^A), where A is the Laplace-Beltrami operator, one assumed that Kt was the kernel of f(t^D), where 
D is any smooth positive second-order elliptic differential operator, acting on C°°(M). It is natural to 
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conjecture then that something analogous holds if D is instead a smooth positive second-order elliptic 
differential operator, between sections of line bundles on M. 

In fact, the following theorem says that this is the case for D = Ag : C°°{'L^) C°°(L'*): 

Theorem 6.1 Let s be an integer. Suppose f G 5(]R+), /(O) = 0. Let Kt he the kernel of f{t^As). 
Then: 

For every R, R' e 50(3), every pair of compact sets J-'r C Ur and J-r' C Ur', and every pair of C°° 
differential operators X (in x) on Ur/ and Y (iny) on Ur, and for every nonnegative integer N, there 
exists c such that for all t > 0, all x G Tri and all y € Tr, we have 

^- («+/+,/) 

\^yK,n-A^,y)\<cj^^^-^, (107) 

where I := degX and J := degY . 

Note The "needlet-type" case, where / has compact support away from the origin, is easiest, and that 
is the only case we will deal with in this article. The general case will be dealt with elsewhere. 
Proof We have 

/(t^A,) = J2 fit^^is)Pis (108) 

Z>|s| 

(strong convergence), so that 

Kt^R'A^, y) - E fi*^^i-)KR',R- (109) 

l>\s\ 

(The interchange of order of integration and summation, needed to derive (jl09p . is easily justified through 
the rapid decay of /, and through ((7T|l . once one recalls ([95]) .) Now, say s > 0. Then by Theorem 15. 6( 



Kt,R,M^,y) = d^nldRy E fit'X,g)f—^K^'ix,y), (110) 

i>\s\ + 

(The interchange of order of differentiation and summation, needed to derive piOp . is easily justified 
through the rapid decay of /, and through ((7T|) in the case s — 0, once one recalls (|95|) in the case s = 0.) 
Now notice that for n G Z, we have 

{l + n){l + l-n)^l{l + l)~j„ (111) 

where 

7n =n{n-l). 

Also note that 

{l + s)l 



= [{I + s){i + 1 ~ sMii + s ~ m + 1 - (s - 1))] ...[{1 + m + 1-1)]. 

{l + s)\ 



il-s)\ 

Accordingly 

^ = [1(1 + 1) - jMI + 1) - Is-i] • ■ • m + 1) - 7i]- (112) 



{i-sy. 

Note also that 

A,, - l{l + 1) - 7_, = 1(1 + 1) - sis + 1). (113) 

To prove the theorem, let us first note that if wc restrict to i e [Tq, Ti], where < Tq < Ti are fixed, then 
the result is trivial (even if we do not assume / has compact support away from the origin). Indeed, since 
d{x, y) < 2n for all x, y, it is enough to show that XYKi^r/m is bounded for (i, x, y) G [To, Ti] x Tr' x Tr. 
This, however, is evident from (I7ip and the rapid decay of /. 
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For other t, we shall assume that the support of / is contained in some interval [a, b] with 6 > a > 0. In 
that case, there is a Ti > such that Kt^R^R' = for all t > Ti (indeed, this will be true if XssT^ > b.) 
Thus we need only prove that there exists a Tq > such that the result holds for all t E (0,To]. In fact, 
choose any To > with 'f^T^ < a/2. For < t < Tq, define 

The function ft is supported in [a, b + s{s + l)ro ], a fixed compact interval. (Note that the denominator 
in (|114p does not vanish for u in this interval; in fact, each factor is at least a/2.) Then, in (|110p . by 
(|112p we may write 

E ^(^'^'^)7n-^^'"(^' = E Mt'>^i)K''{x, y) = t^^K^t^x, y) (115) 

1>S ^ 1>S 

where iCjt] is the kernel of ft{t'^A). Now, it is easy to see that the functions ft {t G (0, Tq]) form a bounded 
subset of 5o(K~''). Thus, by the remarks following (|106p . for every pair of C°° differential operators X' (in 
x) and Y' (in y) on 5^, and for every nonnegative integer N, there exists c such that whenever < t < Tq, 
and x,y E S'^, 

where /' := degX' and J' :— degF'. In fact, this remains true if we insist x G ^r', y G ^r, and allow 
X' (resp. Y') to be smooth differential operators only in a neighborhood of Tr' (resp. J^r), since there 
are smooth differential operators on all of S'^ which agree with these on J^ri (resp. Tr). In the situation 
of (|107p , note that d^^]^ is a smooth differential operator of degree s in x in a neighborhood of Tr/ , and 

— [s] 

dj^y is a smooth differential operator of degree s in y in a neighborhood of J-ri . Thus we have 



\XYKt,w,R{x,y)\^t'' XYd^^,J^yK[t]{x,y) 



< Ct^" —TT: TTTTT^ = C- 



{d{x,y)/t)^ {d(x,y)/t)^' 
as desired. Similar arguments work if s < 0. This completes the proof. 

Remark 6.2 As in the case s — 0, the closed graph theorem implies that the constant c in (|107p can 
be chosen to be < C||/||m, for some M. (Here, for / e 5(]R+), let us define ||/||m = J2i+j<M W^^'^^fWy 
where || || denotes sup norm.) Also, C and M depend only on I,J,N,R' and R.) Note, then, that (jl07p 
remains true if Kt is the kernel of /(t^A^), where now only / e S^'\M.+ ) = {/ G C^(K+) : /(O) = 
and ||/||a/ < oo}. This is evident, since 5o(K+) = nL5^(M+) is dense in 5^-^(K+). Again we will can 
choose c to be < C||/||7\/, where C and M depend only on /, J, N, R' and R. 

Let us now explain what we shall mean by spin wavelets. Say / G 5(M+), / 7^ 0, /(O) ~ 0. Say s is an 
integer, R G 5*0(3), t > 0, and x G Ur. Let us define 

WtxR = X! X! /(^^•^''') sYi„ir{x) sYlm- (117) 
i>ls| m 

Then WtxR e C°°(L''). Also note that 

WtxR.,R' (y) = Kt.R.R' {x, y) (118) 

(if X G Ur, y G Ur', and Kt is the kernel of /(i^As)). This generalizes the case s = of (I104p . Moreover, 
ff F G ^^(L^), then for x e Ur 

Pt,F,.,R ■■= {F, WtxR) = {f{t'A,)F)R{x) = {f3t,F)R{x) (119) 
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where we have set P^. p = (/(t^As)F). We cah fit Fx r ^ ■'P*" wavelet coefficient of F. If F = 
T>i>\s\ Em aim sYim e i^(L''), then 

/^t.F = 5Z .^'/mW. (120) 

Z>|s| m 

SO that 

/3t,F,.,fl = EE fit^^ls)aim sYlmRix). (121) 

i>|s| m 

In cosmological apphcations, one assumes that / is real-valued. In that case, by (|120l and ([50]) - we 
have that 

= [/5t(^)]E; /3t(i^M) = [/3*(F)]m (122) 

for aU F e L^iL""). 

Let us now explain how, in analogy to the case s = 0, one can obtain a nearly tight frame from spin 
wavelets. Let P — P\s\,s be the projection onto the Ti-isi^s, the null space of Ais (in C°°(L*). Then for 
a > 1 sufficiently close to 1, for an appropriate discrete set {xj^k}{j,k)£Zx[i,Nj] on M, certain Rjk with 
Xj^k G C^-Rjfej and certain weights f., the collection of {<t>j k '■— f-j k'^a^ ,Xj k,Rj k} constitutes a nearly 
tight wavelet frame for (/ — P)L^(L'*). (Note that this property is independent of the choice of Rj^k-) In 
that case, we will have that for some C, if e (/ — P)L^(L'*), then, in L^, F ~ C'Ej fc(^j 4'j,k)4'j,k (note 
that this sum is independent of the choice of Rj.k)- We need to make this all precise. 

In order to find suitable Xj^k, we need to subdivide S"^ into a fine grid. To this end, we select co,(5o > 
with the following properties: 

(*) Whenever < r < (5, we can write S"^ as a finite disjoint union of measurable sets £k such that 
the diameter of each £k is less than or equal to r, and such that the measure of each £k is at least cqt^. 

It is easy to see that such r, 6 exist. (We could even require that each set £k be either a spherical 
cap, or a spherical rectangle bounded by latitude and longitude fines. The existence of Cq, 5 as in (*) can 
even be established on general smooth compact oriented Riemannian manifolds; see the discussion before 
Theorem 2.4 of [T^.) 

By using Theorem 16. 1[ one can prove the following result, which generalizes Theorem 2.4 of [19] (when 
the manifold there is the sphere). 

Theorem 6.3 (a) Fix a> 1, and say cq,5q are as in (*) above. Suppose J e iS(M+), and /(O) = 0. Let 

Kt be the kernel of f{t'^As), and let WtxR be as in jllT^ and U18]) above. 

Then there exists a constant Co > {depending only on f,a,CQ and Sq) as follows: 

Let J' C 7j be either finite, or equal to all of Z. Say < b < 1. For each j G J , write as a finite 
disjoint union of measurable sets {Ej^k '■ ^ l£ k < Nj}, where: 

the diameter of each Ej^k is less than or equal to ba\ (123) 

and where: 

for each j with ba^ < So, fi{Ej^k) > co(6a-')^. (124) 
For each j, k, select Xj^k G Ej^k o-nd select any Rj^k G 50(3) with Xj^k G Ur^ k- 

'P.-.fe ^^i'a. eC-(L^). (125) 

For F e L2(L"), set 

jej k 
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Here, if J' = 1^, the series converges unconditionally in L'^iJ^"). 

LetQ^ = X^jGj if J — this series converges strongly onL'^{h'^). Then for all F e L^(L*'), 

I {{Q^ - S^)F, F) I < Cah{F, F) (126) 
(or, equivalently, since Q-^ — 5'^ is self-adjoint, HQ'^ — < Cob.) 

(&) In (a), take J — Z; sef Q = Q^, 5' = . Suppose that the Daubechies condition holds. Let 
P = P\s\.s be the projection onto the 7i|s|,s, the null space of As (in C°°(L'*)j. Then 

{Aa ~ Cob){I -P)<S<{Ba+ Cob){I - P) (127) 

as operators on L'^{Li^). Thus, for any F E (I — P)L^(M), 

(A, - Cofo)||F||2 < ^ ^(i?,- ,)|(F, ,)P < (i?, + Cofo)||F|p, 

SO that, if Aa — Cob > 0, then {/i(£^j fc)-'^/^(p^ ^ is a frame for (/ — P)i^(L*), with frame bounds 
Aa — Cob and Ba + Cob. 

Remarks 1. By ^ and ©, Ba/Aa = 1 + 0(|(a - l)2(log|a - 1|)|); evidently {Ba + Cob)/{Aa - Cob) 
can be made arbitrarily close to Ba/Aa by choosing b sufficiently small. Thus, Theorem 16 .31 gives "nearly 
tight" frames for (/ - P)L^(L'). 

2. The proof of Theorem 16.31 follows along the same lines of the proof for s — 0, which is Theorem 2.4 of 
|19) . It will be given elsewhere [H]. In this article, we will only use Theorem 16.31 in Section 8. 

3. In the statement of Theorem 2.4 (a) of [T9^, was not assumed to be a finite set, but rather, a cofinite 
set (that is, was assumed to be finite), and one had unconditional convergence of the sum for S'-^ F 
and strong convergence of the sum for Q'^. Note, however, that this is equivalent to (a) for finite sets, 
since g-^ = - g■^^ S'^ F = S^F - S'^'^F. 



7 Asymptotic Uncorrelation 

We define a random spin s field G(.) by assuming that there exists a probability space (fi, 3,P) such 
that the map {x,uj) — > G{x,llj) is B{S'^) (8) 5 measurable, 3 denoting a cr-algebra on fl and B{S'^) the 
Borel (T-field of . (To clarify: G{x,ijj) is a.e. in (the fiber of above x), and to say that the map 
is measurable is to say that, for any R g 50(3) the complex-valued function over Ur x ^ obtained from 
the map by trivializing the bundle over Ur is measurable.) As is customary, the dependence of G upon 
uj will be suppressed in the notation. Also, we assume that 



E 



S2 



\C{x)\^dS{x) 



C < oo 



which in particular entails that x — > G{x) belongs to i^(L'*) with probability one (see for instance [40]). 
In cosmological applications it is natural to introduce an isotropy condition: 

Definition 7.1 Let G be a random spin s field. We say that G is isotropic (in law) if for every 
xi, . . . , xn e Uj, the joint probability distribution of Cf{xi), . . . , Gf{xN) is independent of R E 5*0(3). 
Similarly, if G, H are both random spin s fields defined on the same probability space, we say that G, H 
are jointly isotropic (in law) if for every xi, . . . , xjv, yi, . . . , ym € Uj, the joint probability distribution of 
Gf{xi), . . . ,GfixN), Hf^iVi), ■ ■ ■ , Hf {vm) is independent of R e 50(3). 

Note that G is isotropic if and only if, for every R' e 50(3), and for every pi, . . . ,pjv € Ur', the joint prob- 
ability distribution of G^,{pi), . . . , G^,{pn) is independent of i? e 50(3). Indeed, say pj = R'xj, Xj € Uj 
(for j = 1, . . . , TV). One need only then observe that, by ([28]) . G^,{pj) ~ Cf^ (xj) for each j. Similarly, 
G,H are jointly isotropic if and only if, for every R' G 50(3), and for every pi, . . . ,p7v,<Zi, . • ■ ,9a/ € Uu/, 
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the joint probability distribution of G^,{pi), . . . , G^,{pN),G^,{qi), . . . , G^, (qjv) is independent of i? e 
50(3). 

If G is isotropic, then for any R' e 50(3) and any p G Ur', the probabihty distribution of G^,{p) is 
independent of i? G 50(3). Restricting R to the subgroup of 50(3) which fixes p, we see in particular 
that Gr'{p) has the same probability distribution as Gr'r{p) — e^^"^ G ri [p) , where ?A is the angle from 
Pri{p) to Prir{p)- If s 7^ 0, we may take i) with e^^"^ = —1. We then see that the probability distribu- 
tions of ^Gri{p) is the same as that of —^Gri{p), similarly for '<sGr'{p); in particular, both ^Gr'{p) 
and '<sGr'{p) must have expectation 0. Again if s 7^ 0, we may take ip with e^^^ = i, and we see now 
that ^Gri (p) and '^Gri (p) have the same probability distributions. 

Say G is isotropic, and that G^ has spin s spherical harmonic coefficients (for R £ 50(3)). Then, 
evidently, the aj^ all have the same probability distribution, and each has expectation zero. (Here 
= (G^, sYiTn) = {G^,, sYimw) (for any R' e 50(3)) where ( ) is the inner product on 5^.). 

If we take i? to be a rotation about N, we see (from — (G, sYi^ )) ^'^^ from ((60|) that, if m 7^ 0, 
then the probability distribution of aim is the same as that of e^^aim for any 9. In particular, diaim and 
Sflim will have the same probability distributions, if to 7^ 0. 

If G, H are jointly isotropic, af^ = {G^, sYim), bf^ = {H^, sYi„i), then for any I, to, I', to' the covariances 
^("-hnh'm') evidently independent of R. 

In fact, as in the case s = of [Gj, we can say much more about these covariances: 

Theorem 7.2 Suppose G,H are jointly isotropic random spin s fields. Let aim — (G^, syim), and 
him = {H^, sYim)- Then E{aimbi'm') = unless I = V and m = to'. Moreover E{aimbim) does not 
depend on to; we will denote it by Ci^g,h- 

Proof We define Wigner's matrices by 

.Yll^Y.D'n',n{R)sYlm' ■ (128) 

m 

These matrices surely exist and are unitary, since the sYi^ form an orthonormal basis of Tiis whenever 
R E 50(3) and I > \s\. The matrices are independent of s, as one sees at once from (|128p . use of the 
d and d operators, and §7^, ([Ml), (gQl), dSS]). Since .Y^^^' = Em^L'm(^) sY^^, for any R,R', we 
see that in fact the map R ^m'miR) is a unitary representation of 50(3), see also [JS] and [15]. 
This representation is evidently isomorphic to the action of 50(3) on Tiio, which as is well known, is 
irreducible; moreover, the representations for different I are not equivalent. 

Now let af^ — {G^, sYim), and bf^ = {H^, sYim)- For each 1,1' define the covariance matrix M" by 
^mm' = E{aimbi'm') = E{a^bi,m') any R. Noting that af^ = (G, si^;^ ), one easily computes that 
M" = [Z?'(i?)]*Af" D^{R) for any R. The theorem now follows at once from Schur's lemma. 

In cosmological applications, having an "asymptotic uncorrelation" theorem of the following kind is of 
great importance. Generalizing a key result in we now show: 

Theorem 7.3 Let G,H be random, jointly isotropic spin s fields. If x e Ur, let WtxR be the spin wavelet 
of I _?_??] ), where f has compact support away from the origin. Assume that Ci^g,h = 3(0 f^''^ smooth 
function g on the interval (|s|,oo), which satisfies the following condition, for some a > 2: for every 
i G No there exists Ci > such that 

|g«(u)| <c.u-"-\ (129) 
Assume also that, for some c > 0, certain p, t G M, and for all sufficiently large I, 

Ci,G,G > cl-P, and Ci,H,H > cl-\ (130) 
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For X G Ur, let l3c,t,x,R =< G^WtxR >, PH,t.x.R^ ~< H,WtxR >■ Then for any M e Nq, there exists 
Tq > and a constant Cm > 0, such that for all x,y £ Ur, < t < Tq, we have 



l^«K/3c.,..,/3.,...)l<r-^-^^/^(^^^ (131) 



In particular, for fixed x,y, \Cor{l3Q f. ,^ ]:,,(3;[j^^ y p)\ ^0 as t ^ 0+. 

Remark: Note that |Cor(/3j x R^Pty r)\ independent of R. 
Proof By definition, 

iri m Q M \^kf^G,t,x,Rf^H,t.y,R)\ 

\^0AliG,t,x,R.PH,t,y,R)\ = , — / ■ (132) 

^miiG.t,y,R?)^mB^t.yA^) 

Say G, H have spin s spherical harmonic coefficients aim , hm respectively. Note that 

f^G,t,x,R = X! X! fi*^^ls)aim syiniRix), (133) 

i>|s| m 

and similarly for H. Accordingly 

EiPG,t,x,RPH,t,y,R) = E E (1^4) 

l>\s\ m 

= 5] |/P(t2A,,)G,G,Hifl|,fl(a:,y). (135) 

l>\s\ 

To estimate E{\/3q j ^ fll^)i us note that, since / ^ 0, we may choose cq and < a < 6 such that 
l/P > Co on [a2, 6^]. Thus |/p(i2A,,) > cq if a/t < \/A^ < b/t. Choose a', 6' with a < a' < b' < b. Since 
limj^oo = 1> there is a Ti > such that if < t < Ti, \ f\'^{t^\is) > cq if a' /t < I < b'/t. Choose 

< T2 < Ti with (6' - a')/t > 1 if < t < T2. Then, for t sufficiently small, by Theorem [53] and (fT30)) . 
for some C", C > 0, 

E{\f3a,,x,R\')> E cor"^ > ^ ^^2+". (136) 

Similarly, 

EUPhmI') > (137) 

We estimate the numerator in (|132p by using Theorem 16. II In order to do this, we must transform g in 
some simple ways. 

Say s > 0. Reca ll that A;s = (/ - s ){l + s + 1) = l{l + 1) - s{s + 1). Thus, if w > 0, t; = \is if and only 
if / = + + 4[v + s{s + 1)] := Q{v), say. Then Q : (0,oo) — > (0, cx)) bijectively. For v > 0, set 

h{v) — g{Q{v))] then = h{Xis). Moreover, it is easy to see, from (I129p . that for every i e No there 
exists Cj' > such that 

|^M«)I<C>-"/^- (138) 

From p35p . we have 

E{PtxRl3tyR) = E (139) 

l>\s\ 

= rXt,i?,fl(a;,y), (140) 
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where Kt is the kernel of f^t]{t^^s)i where now 

/[.iM^I/pHr-MW^')], (141) 

(smce \f\''{t^\is)h{\i,) = (t^A,,)). Select iX A < B with supp/ C [A,B]. Then supp/[t] C [A,B] 
for all t. By use of l|138p . one sees readily that for any L £ No, there is a C > such that ||/[t]||c^ ^ d 
for all t. Accordingly, by Theorem 16.11 and Remark 16.21 for every M there exists C^^ > with 

— C t"t^^ 

Using this in (|132p together with (|136p . we find the desired result. Similarly for s < 0. 

Remarks 1. Since supp/[t] C [A, B], when deriving the estimate ||/[t]||c^ ^ from (|138p for t small, 
one need only use (I138P for v large (by (|14ip ). Thus, instead of assuming g is smooth on (|s|,cxd), and 
that ()129|) holds there, one may assume only that for some T > 0, g is smooth on (T, oo), and that (|129|1 
holds there. Equivalcntly, we are assuming that g{u) agrees with an ordinary symbol of order —a for 
large u. 

2. By remark (|6.2p . an examination of the proof of the theorem shows that for each M there is an L e No 
such that, instead of assuming that g is smooth and satisfies p29p for all i, we need only assume that 
g eC^ and satisfies ([^5)1 for all z < L. 

3. The hypothesis (|129p . which is used only to estimate the derivatives of h and then the derivatives 
of /[t], can evidently be relaxed - at a cost - to |(7^*^(u)| < CiU~"' whenever i < L, for certain ai. The 
cost in doing this is that one would need to multiply the right side of (|13ip by a factor of for some 
sufficiently large N. It would still evidently follow that for fixed x,y, \Goi{Pq ^ ^ j^, (3jj f y ^ as 

4. If G = H, in applications of the theorem, it is natural to assume a = p = t. 

5. In cosmological applications, where we have (|122p . and where s — 2, one assumes that the polarization 
field is (a single sample of) an isotropic random field F, and that F-e and -Fm are jointly isotropic. 



8 Stochastic Limit Theorems 

In this final section, we give a brief glimpse of the significance of the results of the previous section in 
cosmology. We prove some stochastic limit results, and then indicate why they are important. 
Say e > 0. 

We use a nearly tight frame, as provided by Theorem 16.31 Specifically, in that theorem, we fix a real 
/ supported in [l/a^,a^] for which the Daubechies sum X]jl-oo /^(^^"''") = 1 u > 0, so that 

Aa = Ba = 1- We then choose b sufficiently small that Cob < e. We then produce {(f>j j,} as in Theorem 
16.31 and let j, — fi{Ej^ky/^(pj j,, so that /.} is a nearly tight frame for (/ — P)L^(L''), with frame 
bounds 1 — e and 1 + e. Also, we have (|126p for general finite subsets J^. Of course our choice of ^} 
depends on e; if we want to indicate the dependence on e, we will write 4>j ]^ ^. 

Say G = ^ Aim sYim is an isotropic random spin s field, and let G; = Ci^g,g (notation as in Theorem 
17. 3p . For j e Z, we also let 

7,= Ci{2l + l). (142) 

a-2<a2jAi,<a2 

(To avoid confusion, the summation in p42p . and similar summations in the sequel, are over those I for 
which A/s is defined and for which < a^^ Xis < a^; recall that A/s is defined for I > \s\.) Note that, 
for any I, there at most 3 values of j for which a^^ < o?^Xis < a? , so that 

Xl^j ^ 3 ^ G/(2; + 1) = 3VarG. (143) 

i i>\8\ 
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We let (3jf. = [3jf.^ — {G^cfij f. ,:). We focus on the quadratic statistics 

r.-^El/^.-^r (144) 

fc 

and 

rj= E fia''>^is)\Airn\^- (145) 

/ > I s I , m 

Evidently 

^(Pj) = E /'(a''Ai,)Q(2? + 1) - f{a^'^is)Ci{2l + 1). (146) 

l>\s\ a-2<a2iA,5<a2 

We have: 

Proposition 8.1 E{\fj - f j|) < 67^ for all j. 

Proof For each j e Z, let 

G,= ^imsYi^. (147) 

Then 

i>|s| m 

Note also that (/ — P)Gj = Gj, since no term with I — \s\ can appear in the summation in (|147p (for 
then Xis = 0). By Theorem 16.31 (a) with J' = {j}, we see that, in the notation of that theorem, 

Tj = (S'-^GjjGj), 
f,^{Q^G„G,), 

so that, by Theorem 16.31 (a). 

The proposition now follows at once if we take expectations of both sides. 

In Proposition l8.ll Tj — J2k \(^jke\ depends on our choice of e, and we can even let e << 1 depend on j 
here. By Tchebychev's inequality, the probability distribution of is then a small perturbation of that 
of Tj. We conclude by briefly discussing the latter. 

We will look at random spin s fields G = ^ Aim sYim for which Aim = ^;,-m for all m (so that, in 
particular, Aiq is real). For want of a better word, let us call such a G involutive. Note that, by (|5D|) - 
(I83p . if F is a random spin s field, then both and i^M are involutive. 

We will also be assuming that G is Gaussian, by which we mean that {5RG/(a;) : x G Uj} U {5G/(a;) : 
X e [//} is a Gaussian family. In that case T := {^Aim} U {QAim} is also a Gaussian family. Recall 
that, if TO 7^ 0, and {^Aim} have the same probability distributions, so that they have variances 

El^AimW El^Aiml" = Ci/2. Since we are assuming G is involutive, 'SiAim = {Aim + Ai^-m)/2, 
'<sAim = —iiAim — Ai^-m)/2. It foUows easily from this and from Theorem 17.21 that the elements of 
are pairwise uncorrelated. Since is a Gaussian family they are independent. 
We then have: 
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Proposition 8.2 For all j = 1, 2, 

Var {f,} = E [fj - ETj]" = ^ 2f\a^iXu)iCif{2l + 1) . 

/>|s| 

Say now that for some c, C > and some a > 2, cl~°' < Ci < Cl~°' for I > 0. Then 

^ ^ \i N{0,1) as j -> -oo , 



Var 



—fd denoting as usual convergence in probability law, and N(0, 1) denoting as usual the standard normal 
distribution. 

Proof Write 

r V2f{a^nis)diAi^ for m > 
XhnU) < f{a^^Xis)Aio for m = ; 

[ V2f{a^nis)^Ai„^ for m < 0. 

For each j , {Xim{j)} is a triangular array of independent heteroscedastic unequal variance) Gaussian 
random variables such that 

/> I s| , m 

Var {xLij)} ^ EXUj) ~ {EXLij)V = VHa'' \is){Cir ■ 
Hence we have easily 

Var{f,] = J2 Var{Xf,^ij)}= '^fHa^' hs){Cif 

l>\s\. m- l>\s\. rn 

= Y,2f\a^=\u){Cif{2l + l) . 

l>\s\ 

To establish the Central Limit Theorem, it is then enough to check the validity of the Lindeberg-Levy 
condition ( 10 ) which here takes the simple form 

lim '^^^«"'<"''^'=<"' maxm Var{Xf^) 



^ y maxa-2<a2j;^^^<a2 max^ /'^{a'^J \is){Cl)'^ 

- Ea-<a^.A.<a^ f^a'^ Xl^mY {21 + 1) 

- ^ A^cx, a2." Ea-<a^.A,.<a= /'{a^^ Xls)l ^ ^ ' 



This completes the proof. 
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